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We compute standard model penguin amplitudes in nonleptonic B-decays to light charmless 
mesons using tree amplitude data to fix hadronic parameters. The leading calculation is carried out 
for the a s {mb) penguin contributions from charm quark, up quark, and magnetic penguin loops in 
the NDR and HV renormalization schemes. Power suppressed penguins that are proportional to 
the chiral condensate are also computed using a new factorization formula for these terms, which is 
derived working to all orders in a s (\fmtK). We demonstrate using SCETi that this formula exhibits 
only small perturbative phases and does not have endpoint singularities. Due to our use of data 
to fix hadronic parameters we obtain significantly more accurate predictions for the short-distance 
standard model penguin amplitudes than have been found in the past. Analyzing data in B — > nn, 
B — > Kn, and B — > pp for the penguin amplitudes we find that standard model short-distance 
imaginary parts are an order of magnitude smaller than current measurements, while real parts are 
up to a factor of two smaller with the correct sign. This difference is most likely a consequence of 
long-distance charm contributions or new physics. Constraints on the type of new physics that could 
help explain the data are derived, and used to show that current data favors sizeable long-distance 
strong phases. 



I. INTRODUCTION 



£>-physics experiments have made considerable 
progress in improving our understanding of standard 
model CP violation Q| . Several analyses have fairly small 
theoretical uncertainty and yield precise results, such as 
sin(2/3) from B — > J/^K S or B — > rj'K s type-decays. 
However, for a large number of observables, extracting 
short-distance information depends upon our ability 
to handle QCD effects. Many of these observables are 
sensitive to new physics, and thus considerable effort has 
gone into understanding how to calculate strong decay 
amplitudes with controlled approximations Exam- 
ples of the type of observables are the magnitude and 
relative strong phase of penguin contributions in charm- 
less non-leptonic S-decays, B — ► 7T7t, Kit, pp etc, which 
have significant contributions from loop-dominated 
penguin amplitudes. 

In this paper we classify standard model contributions 
to penguin amplitudes using the SCET factorization the- 
orem for nonleptonic decays from Ref . j3| , and compute a 
missing set of 0(a s (mb)) short-distance perturbative cor- 
rections. (These missing corrections were also recently 
computed in Ref. [ p5| , and we compare results at the 
end.) In principle these corrections have the potential 
of making up for an observed shortfall in explaining the 
penguin amplitude data with leading order strong phases. 
However, we find that these contributions to the ampli- 
tudes are quite small. We also derive a new factoriza- 
tion theorem for "chiraly enhanced" penguin amplitudes, 
which are suppressed by 1 /m& but enhanced by the chiral 
condensate. Our result involves a new generalized form 
factor C? M (z) and a single twist-3 meson distribution 



4>pp(u), and it does not suffer from endpoint divergences. 
We find that these contributions also have small imagi- 
nary contributions. Indeed, all known imaginary short- 
distance corrections to the penguin amplitudes are small, 
roughly an order of magnitude below the experimental 
values in B —> ttt: decays and B — > Kit decays. Expla- 
nations for this discrepancy from long-distance standard 
model contributions are critiqued and weighed against a 
beyond the standard model explanation. 

In the standard model the amplitude for a channel 
B — > M\Mi may be written as 
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and where we use CKM unitarity to 
d or / = s). In this paper T MlAh 
and p A/ i A/ 2 w m be called tree and penguin amplitudes 
respectively. We derive amplitudes for all two-body pseu- 
doscalar and vector modes that do not involve isosinglets 
in the final state. In comparing with experimental pen- 
guin amplitudes extracted from data, we focus on the 
B — ► 7T7T, B — ► Ktt, and B — ► pp channels. 

With the latest data one may extract values for the 
penguin amplitudes in the B — > 7T7t and B — ► pp channels 
using isospin symmetry Isospin implies that P 71 " 71 also 
appears in the 7r°7r° channel, and is absent for 7r°7r~ (up 
to small electroweak penguin terms The same is 

true for PP + p (using the fact that the p's are measured 

to be primarily longitudinal, /£ p ~ 98% and f p L p ~ 
91% 0, and neglecting interference due to the large rho 
width m). To quote experimental values for the penguin 
amplitudes one must pick a phase convention. We take 
T v T and T p 9 to be real and positive, and quote 
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other phases relative to this. For the penguins we will 
quote results for 



p-K-K _ p-n 



pK-n _ p-n 



ppp = ^pp + P f (2) 



etc. In addition, we must also fix the value of the well- 
determined weak phase (3 = 21.2° jl| and the less well 
determined weak phase 7. The latest global CKM fits 
give fi I 



CKMfitter 
/global 



59 g ijo 



UTfittcr 
'global 



64.6° ±4.2°. (3) 



An alternative method to obtain 7 is to use B — > ttt: 
or B — > pp data alone. In principle for £> — > 7T7t this 
is possible using only isospin JlQl , and for B — > pp it 
is possible using isospin, the polarization data, and ne- 
glecting the p width. However the current experimental 
uncertainties need information beyond isospin, such as 
an expansion in Aqcd/™6, necessary to obtain results 
competitive with Eq. (|3|) . An approach with small uncer- 
tainties JTlj ] , which we label the BRS method, augments 
the isospin analysis by using the factorization theorem 
for nonleptonic decays in a specific limited way, namely 
to use Im(T v ° v ° /T" + *~) ~ 0[A/E,0(a s (m b )]. Data on 
Bi(B° 7T+7T-), Br(fl- -> tt-tt ), Br(B° -> 7r°7r°), 
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and C ff + ff - or the analogs for B 



PP 



then used to determine 7. With the latest non-leptonic 
data j|, 12 as summarized by HFAG JjJ, this gives 
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(4) 



where we quote the experimental and theory errors sep- 
arately. With these values of 7, factorization in SCET 
exactly reproduces the observed B Q — ► ir°ir° and B° — > 
p°p° branching fractions. There is also a second solution 
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52.8 
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(5) 



that is, however, disfavored by the additional piece of in- 
formation that the form factor parameter (j > 0. The 
value of 7 from B — > pp in Eq. (Q) has not been quoted 
earlier in the literature, but the analysis is identical to 
that for B — > 7777 in Ref. (llj. Currently the global fit 
values in Eq. (||) and BRS values in Eq. (^) are consis- 
tent with each other at the 1-a level. Suppression of 
can also be studied in a convention 

where Ac is eliminated from Eq. ([j]) however in 

this case charm penguins contribute to the tree ampli- 
tudes which can induce contamination by long-distance 
contributions . We will quote numerical results for the 
penguin amplitudes using 7 = 59° and 7 = 74° to give 
some indication of the spread of possible values. 

With the latest B -> tttt and B —> pp data §, [jj), 
the isospin formula quoted below in Eq. ( p"5| ) gives the 



penguin amplitudes for 7 = 59°, 

10 3 P™ = (1.77 ± 0.73) - i(2.91 ± 0.58) . 
10 3 Ppp = (-2.91 ± 2.63) - i(0.78 ± 1.82) 

while for 7 = 74° we find 

10 3 P 7 ™ = (4.41 ± 0.61) - i(2.91 ± 0.58) , 
10 3 P pp = (3.81 ± 2.34) - i(0.78 ± 1.82) . 



(G) 



(7) 



Here for convenience we pulled out a prefactor to quote 
a dimcnsionless penguin amplitude P, using 



P M * M * _ G F m 2 B 
(lGcV) - V2 



P 



M1M2 



(8) 



Note that for fixed CP-asymmetries C and 5, the ex- 
tracted real part of P 7r7r and _P PP depends fairly strongly 
on the value of 7, but the imaginary part is completely 
independent of the choice for 7. (This is demonstrated 
explicitly in Eq. (|l5|) below.) Neither result depends on 
the error in \ V u b\- The experimental errors here have de- 
creased noticeably from early penguin extractions [fl4|| . 
The challenge for standard model predictions is to repro- 
duce or rule out the values in Eqs. (|6|-[7]). 

The extraction of both the real and imaginary part 
of penguin amplitudes in the Kir system currently re- 
quires further theoretical input. In B — > Kix decays the 
penguin amplitudes dominate the tree amplitudes due to 
CKM suppression, making a precise comparison of their 
values even more interesting. Both types of amplitudes 
are important in CP asymmetries. Using a A/m^ expan- 
sion, the tree amplitude for B° — > A' _ 7r + at LO depend 
only on hadronic parameters C B7r and (j n that are fully 
determined by the tree amplitudes in the B — > ttt: chan- 
nels, plus the cf>K twist- 2 distribution function fl5|| , This 
allows the phase of the penguin amplitude P K+7T to be 
extracted from the data using only factorization for the 
tree amplitudes (which we will refer to by adding a sub- 
script TF). The tree amplitudes are reliable since a proof 
of factorization to all orders in a B was given in Ref. j^. |T(i| . 
extending the original proposal and one-loop analysis in 
Ref. JTHj . Although factorization has also been demon- 
strated for light-quark penguin loops (it, d, s), a complete 
analysis for charm-loops is still lacking. Using the phase 
convention where T K 11 is real we find for 7 = 59° 



10 3 P£ F 



while for 7 = 74° 



(4.87 ± 0.39) -i(2. 22 ±0.77) 
-(4.22 ± 0.36) - i(2.22 ± 0.34) 



10 3 P, 



K~ 
TF 



(4.73 ±0.36) 
-(4.41 ±0.34) 



i(2.16±0.73) 
i(2.16±0.68) 



(9) 



(10) 



The only B — > Ktt data used here was Br{K~ tt + ) and 
Acp(K~tt + ), and there are two solutions for each 7. Al- 
ternatives to the above analysis extract the Ktt penguin 
amplitudes using a SU(3) based analysis with the tttt 
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data |ll| or by a global SU(3) based fit |2(J, and these 
yield similar conclusions for the size of the penguin am- 
plitudes. Again the data gives i^7r-penguin amplitudes 
with large imaginary components which require explana- 
tion in the standard model. 

To determine penguin amplitudes for charmless B- 
decays in the standard model it is convenient to orga- 
nize the relevant mass scales as an expansion in A/m b 
and A/m c This can be done from first principles 

using the effective field theory SCET . In this expan- 
sion certain contributions to these amplitudes factorize 



allowing them to be parameterized by well defined uni- 
versal hadronic matrix elements. Since we are interested 
in the standard model prediction, wc also organize the 
amplitude according to large (61,2,85) and small (6*3-10) 
Wilson coefficients. To explain which terms will be com- 
puted in this paper we schematically give a result Po for 
channels B — > MM' with pseudoscalars, MM' = PP, 
with pseudoscalars and vectors MM' = PVq, and with 
two longitudinal vectors MM' = VqVq. For complete- 
ness we also quote an analogous result Pt for transverse 
polarizations MM' = V T V T : 



P> 



(11) 



P T 



65,6 + 



a s (mb)Cx,2,Sg\ \HW ^BMj.M' , MM' t BM ^M'l ! (c Ai \ as ( mb ) Cl - 2 ^9 \ ^M_ qBMiM' 

BMM' , a s( rn 'b) ( n f_xMxM' 1 n t_x+J$MiM'\ , r < Us ( m b)^M t A M xM' 



MM' i-BM ±M' , MM' r BM ,M' 
S> H>pp T — — t, j <P„ P 

m b FF m b yy 



+ C h2 a s (2m c )vAf 5 



A <~ p 



m b 



^ppi 



(12) 



+ 63,4 + 



a s (m b )Ci t2 ,8g 



~ «s(m b )Ci,2,8 g ^ 1 >BV,K' 

°5,6H j W_L ( f , _L 

7T / TOfc 



BW , ^ a s (^b) 



Ci :2 a s (2m c >^f/ V + 65,6 



.4 



svv' 

(2ann — x) 



When coefficients 63,4,5,6 appear here we leave implicit 
the fact that the electroweak penguin coefficients 67,8,9,10 
can also appear. 

The terms on the first line of Eq. ( pi] ) for Po are lead- 
ing in the A/nib. c expansion. Working to all orders in 
a s (/j,i) at the intermediate scale /i; ~ y/n%bA, they in- 
volve a so-called soft form factor ( BM and a hard form 
factor Cj M - (Though it should be emphasized these 
names are somewhat misleading, since both form fac- 
tors involve hard-collinear exchange, and thus the same 
length scales.) The terms on the second line of Po are 
the "chiraly-enhanced" power corrections suppressed by 
HM'/mb- Here fiM' ~ O(A) is a ratio of the squared 
meson mass to a sum of quark masses, and is important 
for nonlcptonic decays because it is numerically enhanced 
Hm 1 — 2 GcV |l7j . In Eq. (|ll]) we display a new result 
that we will derive below, namely that to all orders in a s 
the chiraly enhanced terms in Po arc determined by one 
new form factor, C^ M an d one twist-3 distribution 
for longitudinal polarizations/pseudoscalars. 

On the third line of Eq. ([ll]) for Po we have a term 
Af s MM . This is the so-called "char ming penguin" due 
to long-distance charm loop effects j22l H| > whose leading 
contribution is expected to come from the charm thresh- 
old region H [24]]. It is parametrically down by a power 
of the non-rclativistic velocity v relative to the lcading- 



M 
PP 



power result. 

The remaining terms on the third line of Eq. (Jll]) are 
due to annihilation and are suppressed by one or more 
powers of A/m b . Terms at order l/m b contribute only to 
longitudinal polarizations and include: fB<fi M <t> M which 
was studied in Refs. p 



M 27 



and the remaining 
leading annihilation amplitude fB4>~B4 >3M 4 >M > which was 
computed recently in Ref. pq] . The [i m fB(f>pp(l) AI ' terms 
are chiraly enhanced annihilation studied in Refs. |27], 
|29|] . At lowest order in the a s expansion the annihilation 
terms shown in Eq. ([ll]) are real |2{| (using the zero-bin 
procedure [p0|). A nonperturbative complex annihila- 
tion amplitude involving soft exchange occurs at order 
a 2 s (i±i)/m b . 

The terms shown in the Pp amplitude follow a similar 
notation to that for Po, but there are no analogs of the 
LO terms on the first line for Po |3l]] . The term shown on 
the first line of Pt comes from j-p conversion |3^] , and is 
suppressed by a em , but enhanced by m b /A. Pt does not 
contain chiraly enhanced terms, so we show on the second 
line the terms from the analogous operators, which gen- 
erate three form factors £f v , Cj± , an d Cfl- In Ref. @ 
analogs of the terms on the second line of Eq. (jlj) were 
computed with an expansion in a s (Hi), and the penguin 
annihilation term A B ^^ l _ was also treated. Following 
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Rcf. |^4j wc included a term A B y y in the amplitude to 
produce transverse vector mesons. This is the most con- 
servative approach given that so little is known about the 
factorization properties of A cS . 

We remark that the second and third lines of Eqs. ( p4|) 
and do not contain the complete set of 1/ m& or l/m| 
power corrections, but rather a collection of terms that 
are believed to be important due to numerical enhance- 
ment. 1 The additional uncertainty from missing l/m,b 
corrections will be taken into account in our final error 
estimates. Also the C34 and C$fi terms always come to- 
gether with an ot s C\ t 2 term, which is the matrix element 
responsible for canceling the largest scheme dependence 
in these coefficients at next-to-leading-log (NLL) order. 
Since tree level 63- 10 terms could compete numerically 
with a s Ci.2,8g terms, we will require both to be included 
in what we call our leading order penguin amplitude. 
In this paper we will neglect a. s (rn,b)C3^ terms relative 
to a s (rrib)Ci,2,8g since numerically they are 6-30 times 
smaller. This is the same strategy that was adopted for 
the "NNLL computations" in B -> X s £+r §|. 

Eq. ( p"l| ) is schematic because we have not yet displayed 
the precise coefficients in front of each term. The coef- 
ficients for the C, BM terms on the first and second lines 
were computed in Ref. jl7|]. To obtain Eq. (11), we ex- 
panded in A/m and a s (rrib), but avoided the additional 
uncertainties from expanding in a s (ni), where the inter- 
mediate scale /Uj ~ V EA ~ 1.3 GeV. This is made possi- 
ble by the fact that the form-factor parameters (f M and 
£-bm are umversa l hadronic parameters when we distin- 
guish m 2 > EA > A 2 j|. Without expanding in a fl (//*), 
the third term on line 1 was obtained in Ref. §. The 
third term on line 1 and third and non-a s (mb) terms on 
line 2 were computed at leading order in an expansion in 
a s (fJ.i) in Ref. Q. 

A main goal of this paper is the computation of the 
a s Ci a # g C, BM <j) M ' terms on line 1 of Eq. (|ll|). Note that 
this term is a leading-order contribution to the penguin 
amplitudes due to the hierarchy in the Wilson coeffi- 
cients. Our other main goal is to derive the factoriza- 
tion theorem for the terms on the second lines of Pq , and 
compute the corresponding Wilson coefficients. We de- 
rive the factorization theorem working to all orders in a s , 
and perform tree level matching for all contributions. We 
also compute the a s ( BM terms, however the perturbative 
matching computation that determines the a s (mb)Cj AI 
and a s (nib)( BM terms on line 2 will not be considered 
here 

We would like to compute the magnitude and phase 
of P for the 7r + 7r~ , if + 7r~, and p + p~ channels from the 



1 Also note that to our knowledge it has not been demonstrated 
that the chiraly enhanced terms shown in Eq. (hlh give the com- 
plete set of such enhanced contributions at this order in the power 
expansion. In particular it remains undetermined whether time- 
ordered product terms appearing at this order are or are not 
chiraly enhanced. 



terms in Eq. ( pr| ) and compare with Eqs. (|§-[l0|). The 
phase here is that of P/T. A schematic expression for T 
is given by 

f „ (c 1 , 2 + as{mb)Cl ' 2 )c BM 



+ (C 1 , 2 + as{mb)Cl ' 2 )cf M M ' 



a s (m b ) rpBMM 1 



- (ann) 



(13) 



where we suppress terms from Wilson coefficients 
and power suppressed terms other than annihilation. 
For T- IT + 7r - , T p + p - , and Tk+k- the leading term in T 
is real and numerically dominant. Therefore in P/T 
the C\^ol s corrections from P arc leading, while those 
from T are higher order in a s . In fact it would be in- 
consistent to keep the Ci^ag corrections in T without 
keeping the Cs^otg terms in P, because both of these 
terms carry ^-dependence that cancels that in ^ BM and 
Cj M - Thus to compute the imaginary part of P/T 
at the order we are working, we can take J ,7r+7r and 
f K+ ^ to be real. Numerically Refs. §4], || found 
that including one and some two-loop corrections gives 
10 3 f ff7r = (31+1 - i0.07± 3 7 ) + fF^ n) . Thus, the imagi- 
nary part is significantly smaller than the real part. Us- 
ing the zero-bin procedure, the annihilation contributions 
^(ann) 2 are a ^ s0 rea ^ a ^ leading order in an expansion of 
a s at the hard and intermediate scales (2^, [2{| . 

Besides Cj M ( z ) an d C BM > the other hadronic parame- 
ters in Eq. ( |ll| ) include the twist- 2 distribution <f> AI and 
twist-3 chiraly enhanced distribution <\>^ v defined below 
in Eq. (p6|). In this paper we adopt the point of view 
that Q B and the normalization 



r BM _ 
SJ = 



dz{f M (z) 



(14) 



should be fixed using other data (tree amplitudes and/or 
form factors) and then used to make predictions for the 
penguin amplitudes We will see that, relative to 
adopting models for all the hadronic parameters, fit- 
ting to tree amplitudes removes the dominant hadronic 
uncertainty in the computation of the short-distance 

penguin amplitudes. In a generic new physics model, 
poxpt _ psM + pBSM^ go tQ tcst the data for new p h ys j cs 

we must have control over P SM . 

The plan for the paper is as follows. In section || 
we give formulas for determining the penguin ampli- 
tudes and the soft and hard form factor parameters from 
the data. Section |n reviews the leading factorization 
formula, and section IV] discusses the endpoint behav- 
ior of Cj M (z). In section [v| we give a summary of all 
0(a s (mb)) one-loop hard coefficients at LO, and then in 
sections |\T and |VII| provide more details of their calcula- 



tion in the NDR and HV schemes respectively A factor- 
ization theorem for chiraly enhanced penguins is derived 
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in section VII] working to all orders in a s at the interme- 
diate scale \/m,bA. In section IX we discuss long-distance 
charm contributions. Penguin annihilation contributions 
are reviewed in section KL Our analysis strategy is out- 
lined in section 
in section 



XII. 



XI and input parameters are summarized 



Our numerical a nalysi s for standard model pengui ns is 
taken up in section XIII . This is followed by section XIV 
where we derive constraints on the effect of new physics 
contributions, and discuss what is needed to shift the 
penguin amplitudes closer to the data. Fur ther discus- 
sion and conclusions are given in section XV . Several cal- 
culational details are relegated to appendices. 



II. DETERMINING PENGUIN AMPLITUDES 
AND THE C BM AND £f M FORM FACTORS WITH 
NONLEPTONIC DATA 



The B — ► 7T7T data can be used to extract the pen- 
guin amplitude P" w and the tree amplitudes T" T 
and r° I_ , including the strong phase in p* + *~ j T ^^ _ 
Solving equations in Ref. Jll| with our phase convention 
the penguin amplitude is 



Re(P^) = N ^ 



(1 GeV) 



(i™) 2 [sin2/3- sin(2/3+2 7 )] - R c ( sin2/3 + 



V cb V* d \{2t™^)cos2f3 



Im(P™) = -^gZ 
v ; 1 GeV 



(lGeV) |y cb y c *J(2i-sin 7 ) 



(15) 



where the parameters on the right-hand-side are determined by nonleptonic data: 



64tt Bt(B- 



r B - 



S1117 



1/2 



l + Bnn cos 2/3 + SV+tt- sin 2/3 



c 2Bt(B- -> 7r 7r-)r B o ' 



Bw — J 1 C 7r+7r _ ^tt+tt- 



(16) 



Eqs. ( |15| ) and (16) also determine the penguin for longi- 
tudinal B — > pp decays, by simply taking all superscripts 
and subscripts 7r — > p, and were used to determine the 
numbers quoted in Eqs. (|^-0) with \V c b\ = 0.0417 and 
\V c d\ = 0.227. The Im(P) is mainly sensitive to the di- 
rect CP-asymmetry. Since (t™sin7) and (i pp sin7) do 
not explicitly depend on the weak phase 7, the same is 
true for the values extracted for Im(P*" 7r ) and Im(P pp ) 
(demonstrating the statement we made in the introduc- 
tion). The amplitude parameter = |T 7r+7r_ |/|T 7r ° 7r+ 1 
gives information about the size of the color-suppressed 
tree amplitude. The results in Eq. ( |l5| ) and ( |l6| ) are based 
on isospin symmetry and neglect small electroweak pen- 
guin contributions in B~ — > ir°ir~~ . This analysis leaves a 
± sign ambiguity in front of the B nn dependence in , 
which we resolved in Eq. Jl6| ) by taking the "+" solution. 
This solution is preferred by the standard model and 
rigorous power counting for the QCD amplitudes. The 
other experimentally allowed solution, B V1V — * —B wl7 , 
has very large penguin amplitudes, 10 3 |P lr7r | ~ 11 and 
10 3 |P p J ~ 24, which are extremely difficult to accommo- 



date in the standard model. 

An important question for the phenomenology of 
charmless nonleptonic B-decays is the relative size of the 
form factors ( BM and (f M (defined below in Eq. @). 
Here we follow Ref. pi and organize the expansion ac- 
cording to Cj M ~ C B which is a natural power counting 
when factorization is not used at the intermediate scale, 
and is also the scaling used in the KLS approach [ p5| . 
This counting is supported by the B — > rnr data and 
factorization with the zero-bin procedure of Ref. J3TJ| 
(which implies that both of these form factors have an 
a s (pi) at leading order in the expansion at the inter- 
mediate scale). In the BBNS approach JijJ, a hierarchy 
Cj M ~ a s(Hi)( BM is adopted. This changes the order of 
terms in the perturbative expansion which wc discussed 
in Eq. (fLlT) by making certain terms higher order in a s . 



At leading order the B — > mr factorization theorem 
for the tree amplitudes can be used to extract the nor- 
malization of the soft and hard form factor parameters. 
Expressed in terms of obscrvables the result is p6| 
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\V V *\(C B * + Cf*) = 

\V ub \{x-%^ = 



U\v ud \ 



h\Vud\ 



{C 1+ C 2 )tl 



Co — Cft 



3Ci + C 2 



- (C 1 + C 2 )(C i -C 3 
-C 3 +3C 4 -4(Ci + C 2 )^ 



l + 0(a s (m b ), 



Cf 



cl 



(C 1 + C 2 )(C i ^ 

I 



■C 3 ) 



(17) 
(18) 



where E is the pion energy in the CM frame of the B. 
On the left-hand-side of Eq. ( |l7| ) we have the semilcp- 
tonic B — > 7r^j> form factor at q 2 ~ 0, which is given 
by /f^fO) = C B7r + Cf 7 "- The expression for {x" 1 )^ 
in Eq. ([18]) follows in a straightforward manner from re- 
sults in Ref. j|, but to our knowledge has not been pre- 
sented in this simple closed form in the literature. For 
the hadronic parameter (a;" 1 )^ = j^dx (j> 7r {x)/x 1 a fit of 
the 7-7T form factor to 7*7 — > ir° data gives p7| 

(x- 1 )^ = 2.9 ±0.4. (19) 

Using the latest experimental data jjj, 7 = 67 ± 10°, 
\V ud \ = 0.9738, and LL Wilson coefficients, Eq. (Il7|) gives 



/r(o) 



0.182 ±0.011 ±0.036 



ff Cj = 0.262 ± 0.052 ± 0.052 



4.2xl0~ 3 

I Kb I 
4.2 xlO" 3 

1Kb] 



(2°) 



We emphasize that these results do not rely on Eq. (19). 
In Eq. (|20| ) the first errors are experimental and the sec- 
ond theoretical. Theoretical errors are computed as a 
generic 20% error on the central value, both here and 
below in Eq. (|3|). Setting \V ub \ = 4.2 x 10~ 3 and using 
(x' 1 )^ = 2.9 in Eq. (poh gives 



( Br! = 0.092 ± 0.027, 



0.090 ±0.018. 



(21) 



These values favor £j ?7r ~ C Sjr . There is a sizeable corre- 
lation in their quoted errors, and to take this correlation 
into account in our numerical analysis we will do Gaus- 
sian scans over the range of experimental errors quoted 
111 Eq. @. A 12% error in |Kb| is also included in our 
final results. 

The same results, Eqs. ( |l5||l8| ) apply for B —> pp 
for longitudinal p's, where now one uses N p o p - and de- 
termines t p c p from the B — ► pp branching ratios and 
CP-asymmctrics S p + p - and C p + p -. Here the analog of 



/+ (0) is the longitudinal B 



A|,(0) = 



A a (0) 



tv form factor at q 2 
(m B + m v ) 







2mi, 



Ax(0) 



2my(m B +TOy) 

Taking 7 = 67 ± 10° the nonleptonic data gives 

4.2 xl0~ 3 



(22) 



Af'(0) 



0.261 ± 0.022 ±0.052 



1Kb 



" 1 ) p C,f p = (o.06±0.11 



4.2x10" 
I Kb I 



(23) 



Again the first errors are experimental and the second 
theoretical. Due to the large uncertainty in the central 
value of (.t~ 1 ) p £j P a 20% theoretical uncertainty would 
not be noticeable. 

It is interesting to make a comparison of the tt and 
parameters, which from the results in Eq. (|2(]) and (|22 
give 



Bp 



R 



0.23 ±0.42. 



(24) 



This large 0.42 experimental error is induced by our cur- 
rent knowledge of the B — > pp CP-asymmetries together 
with greater sensitivity to 7 and t c . Now in SCETn one 
can derive a factorization theorem for Cj M (discussed in 
Eq. ([l6| ) below) that implies that 



B>pTT 



l + Ola,(jJn), 



(25) 



Here the theoretical errors should be increased to ~ 35% 
to account for the additional expansion in a s at the in- 
termediate scale //.,-. Using f p ~ 1.6 we find that the 
nonleptonic data plus factorization at the intermediate 
scale currently implies 



0.38 ±0.35 ±0.13. 



(26) 



where the first error is experimental and the second is 
the 35% theoretical error. 

A result for this ratio can also be obtained from the 
factorization theorem for the color-suppressed decays 
B° -► D°p° and B° -► D tt derived in Eq.(69) of 
Ref. j$8|, which gives 



'Br(S° -► L>°p°) 
Br(5° L>°7r°) 



0.62 ± 0.24. (27) 



It is quite interesting that this ratio is found to be less 
than unity, and that the results extracted in Eq. ( pjjj ) 
from charmless decays, and in Eq. ( p7j ) from charmed 
final states, agree within errors. The significant range al- 
lowed by the errors can be reduced by noting that there 
is a rigorous lower bound on the inverse moment for pos- 
itive definite <j) M {x), 



{x )m 



dx fl^l> [ dx<t> M {x) = 1. (28) 
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For M = 7r and M = pwe have M (x) = 4> M (\ - x) 
from isospin and charge conjugation, and this bound can 
be strengthened: 



sented in section XII. The basis of operators is 



x 



,1/2 

> 4 / dx<p M (x) 
Jo 



dx ( — 

o v ^ 

2. 



1 



1-a; 



(29) 



With the mean value in Eq. (|T^) this bound is close to 
the central value in Eq. (27). 



For our analysis we take into account the bound and 
the data in Eqs. (|2(]) and (|27|), and hence use a model 
for 4> p {x) that is constrained such {x~ 1 ) p = 2.2^Q2 which 

-3 



gives (x l ) p /{x x ) 
we then find using Eq. (|23j) that 



0.76. Taking \V ub \ = 4.2 x 10 



C Bp = 0.234 ± 0.065, Cj P = 0-027 ± 0.049 . (30) 

Again we scan over the range of experimental errors 
in Eq. (23) to take into account the sizeable correla- 
tions. Note that £j p is sensitive to the values of the 
B — > pp branching ratios, which dominate the error 
and favor a smaller color suppressed amplitude than in 
B — > 7T7T. The central values in B — ► pp are consistent 
with £j P < a s (pi)C, Bp . However, within errors the scal- 
ing that we adopt, (j p ~ ( Bp , is also consistent (given 
that ~ still means that the factors can differ numerically 
by a factor like 2). As the experimental uncertainties 
on the nonleptonic decays decrease, we expect the com- 
bined analysis of B — *• tttt, pp introduced in this section 
to play an important role in furthering our knowledge of 
hadronic parameters appearing in the factorization the- 
orem for charmless nonleptonic decays. 



III. FACTORIZATION AT LEADING POWER 

In this section we review the SCET factorization anal- 
ysis at leading order from || to setup our notation. The 
decays B — > M\Mi arc mediated in full QCD by the 
weak AB = 1 Hamiltonian, which for AS = reads 



G, 



10,77,83 



* p—u.c i—3 



(31) 



where the CKM factor is A 
at LL order 



VpbV*f with f = d and 



Ci_ 10 (m 6 ) = {1.107 , -.249 , .011 , -.026 , .008 , -.031 , 

-4 n -7wm-3 l.gxlO" 3 }. (32) 



4.2 x 10" 4 , 4.2 x 10" 4 , -9.7 x 10 -3 . 



The coefficients in Eq. ( |3~l| ) are known at NLL order |?9| , 
and the values we used for our main analysis are pre- 



0{ 

o 5 
o 7 

Os 
Og 

o 7l 



(pb)v-A(dp)v-A, 

(Ppba) V-A {d a Pf3 ) V-A , 

(db) V -A(qq)v-A , 
(df3b a )v-A(q a qi3)v-A , 
(db)v-A{qq)v+A , 
(dpb a ) v-a {q a qp ) v+a , 

^Y(db) V - A (qq)v+A, 
-^-(dpba) v-a (q a qp ) V+A , 
^Y(db) V - A (qq)v-A, 
-^-(df}b a ) v-a (q a qp ) v-a , 



(33) 



gm b 

em b - 
'4^2 



da^G^Pnb, 



^da^P R b. 



Here a, (3 are color indices and e q are electric charges and 
the q are summed over the light quarks, q = u,d, s, c, b. 
The AS = 1 Hw is obtained by replacing (/ = d) — > 
(/ = s) in Eqs. ( |3l| , |33| ). The numerical dominance of 
Ci, C2, Cs g will allow us to simplify the calculation since 
we need only include the effects of Oi,2,8g at one-loop. 
Perturbative corrections due to the other operators are 
numerically tiny. Our sign for g is such that the QCD 
fermion Feynman rule is igT A "f^ . 

The matching onto SCET occurs in two stages. First 
one matches onto SCETi by integrating out fluctuations 
at the scale m b . One then matches onto SCETn at the 
scale V Am b . For the LO factorization theorem for non- 
leptonic B-decays this second step of matching can not 
lead to strong phases, as discussed in ||, and for ( BM is 
known at one- loop order 40, pd|. In this paper we wish 



to complete the 0(a s ) matching for the first stage. For 
tree amplitudes the corresponding computation was car- 
ried out in Ref. JO. Here we consider the result for the 
penguin amplitudes. In particular we present the short- 
distance up and charm loop contributions in two different 
regulation schemes for 75, as well as corresponding con- 
tributions from the magnetic gluon operator. 

At the scale p ~ m b the Hamiltonian in Eq. (|3l]) is 
matched onto operators in SCET. Due to the nature of 
the matching onto SCETn the first two orders of the 
power expansion of the Hamiltonian in SCETi are needed 
to determine the leading-order amplitudes 

" ^ n, n ^ i 

+ J2 [ [^UbP^Mfi^) + ■■■)■ (34) 
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The operators for the AS = transitions are a. 16 



and it was convenient to define 



and 



Qi 



(0) 



4P L i 



Q2d,3d = [dn^fiPLbv] [un^rf-PL^Un^, 
Qfdfid = [dn^^Phbv} [qn^^P^RQn,^ 



(35) 



(39) 



1 &,<^ P *G 



U>3 

n,ui 3 



i = 1,2,4,5 
i = 3, 6 



Qld = — ig^ u P L b v ] [dn^PLUn^s] , 
TOfc 
r) 

771 b 

TYl}y 
<ry 

Qsd,6d = [dn,ui ig$n,u; 4 PLbv] [qn,w 2 iPL,Rqn,u 3 ] 

Till) 

-2 



(36) 



-2 



Qsd = — -[oV-i igBn.^PLK] [dn, U2 $^PRqh,U: 

TTl}) 

Here P £ = (1 - 7 5 )/2 and P R = (1 + 7 5 )/2. At lowest 
order, Qjdsd §^ ve a vanishing contribution to the rates. 

Q&ed wn l n °t be relevant in our analysis since we will 
not be considering isosinglet final states. At tree level 
the matching onto Q^Jed was done in Ref. p3| . From 

Eqs. (|35| , ^6| ) the AS = 1 operators are obtained by 
swapping d — > s. The "quark" fields in Eqs. (|35|,[56|) with 
subscripts n and n are products of collinear quark fields 
and Wilson lines with large momenta Wj. In particular 
we have defined 



Un )W = [^Vn^W-fi-pt)], 



(37) 



= 7~T [Wj[m.D c ,„,iX)S ) jW n <y(w-P t )] 

where £^ creates a 77-collinear up quark or annihilates 
an antiquark. The b v field is the standard HQET field. 
For a complete basis we also need operators with octet 
bilinears, T A ® T A , but their matrix elements vanish at 
LO. 

The leading-order amplitude is generated by time- 
ordered products of both the operators Q' ' and Q' 1 ) 
with insertions of a subleading Lagrangians G3, pjj . T- 
products with can be factorized as T[Q™ and con- 
tribute to terms with C, BIiI , while T-products with Q^> 
can be written as TjQ™ and contribute to terms with 
( BM . Here 



Tl [Q^} = jd±ydyT[Qf\0) i£W q (y) 

+ /d^dV T[Q<°>(0) iC^fo) i£«(y')}] 
+ fd*yT[Ql°\0),ii:£:?(y)], 

T 2 [Q«] = Jd 4 y T[Qf 5 (0), «41g(^)] ' ( 38 ) 



The Lagrangians in Eq. ( |38| ) can be found in Ref. J46[ . 
Note that only the n-collinear fields appear in the T- 
products T\ and T 2 , which explains why the same T\ 2 
appear for heavy-to-light form factors at large meson en- 
ergies The form factors simply do not have the ex- 
tra Qf. In addition we have opcrators/T-products whose 
matrix elements give A cc . We refer to section IX below 
for further discussion of these contributions. 

In this paper we use factorization at the scale rrib, 
where the hadronic parameters are defined by matrix ele- 
ments of T\ and T 2 and the n-collinear operator, namely 

(M n |Ti [qtM \B) = Cl M 5 Ul m B ( BM , (40) 
(M n \T 2 [q L null i9VLM \B) = -Cl M 8^J^^ M {z) , 

(Mn\q£,JqLM = \ C ™J^ hi M«). 
(Mn\q£jqL 3 \0) = \ C^J^ 3 f M 4>m(u), 



where z = u>i /ttib, u = 0^2/777 _b and we have made the 
momentum-conserving 5- functions explicit, 6^ = <5(wi — 
m B ). <L lW4 = K^i + u i - m s), and 6%% = S(uj 2 - w 3 - 
tob). As pictured in Fig. u and 1 — u are momentum 
fractions for the quark and antiquark 77-collinear fields, 
and z and 1 — z are the momentum fractions carried by 
the n-collinear quark and gluon field in Qif '. Finally, 
Cf M and Cf 1 are Clebsch-Gordan coefficients. We fix 
the following sign convention for the states 



7r + = +ud, 7T° = —=(dd— uu), tt — —du 1 
v2 

K° = sd, K~ = -su, K+ = us, K° = ds, 



B" = bd, B~ 



-bu, B + = ub, B° 



db. 



(41) 



and take vector meson states to have a negative sign rela- 
tive to the corresponding pseudoscalar mesons. The over 
all phase convention is fixed so that the Clebsch-Gordan 

and C p d = = +1. One can then compute that 

Cf^ 71 ' = +1 and Cf rU = etc. Note that the signs 

take into account whether the operators have left or right- 
handed quarks. Putting the pieces together gives the 
leading order factorization theorem which integrates out 



9 



hard ~ m 2 , fluctuations 
A LO = -i (MiM 2 \H w \B) 
GFm 2 B 



(42) 



V2 
+ ( 



-f. 



Mi 



dudz Tu(u, z)Cj Ah {z)(j) Ml (it) 



BAI; 



duT lc (u)(j> Ml (u) 



a 



2) 



Here the hard coefficients Tiq and T\j depend on channel 
specific linear combinations of the matching coefficients 



ri C («)=c. 



BM 2 r M t df) 



pBM-2 pM ± 



(«)■ 



qBM-2 

fh 



r n 'i M) 



IBM; 



(«) 



7 Ml 



Tu(u,z) 



C» M *Cf^\u,z) 



iBM; 



C^ u bi f \u,z) 



Cf^C^A f) (u,z) 



r BM 2 r M ± . ' 



(u,z). 
(43) 

Results for these T's in different decay channels can be 
read off of Table I in Ref. ||. Power counting implies 
(™ ~ Cf M ~ (A/m 6 ) 3 / 2 while </> M - 1. Here the non- 
perturbative parameters ( BM , (f M (z), and 4> M (u), all 
occur in the B — > M scmilcptonic and rare form fac- 
tors. For a model independent analysis they need to be 
determined from data. Note that in the leading order fac- 
torization theorem all terms involve a form factor times 
a meson distribution function. 

Taking the terms proportional to a£^ from Eq. (E3) 
generates the penguin amplitude terms on the first line 
of Eq. (pi]). Using Eq. (|13) still requires matching the 

full theory Oj's onto the q[°/ ^ to determine the Wilson 



( f) ( f) 
coefficients c\ and b\ . 

f = d, s we have 



For the coefficients of with 



df) 



,(/) 



,(/) 



df) 

c 4 



9i 

Ci 

' N, 



9i 

Cio 



( n _i_ ^ 3 ^ w ^9 

V U + W r T 2N r 



Ac\ 



(/) 



Ac 



(/) 



(/), 



,(/) 

4 



and for the Q^ 1 ' we have the coefficients 

-A 



df) 



df) 



df) 



<h + [I- 
C 2 + (l 



uJ N, 
u 



(44) 



(45) 



df) 



df) 



df) 



1\3C 9 
u) 2N C 
Ci 
N c 
1\ 3C 10 



Ab\ 



(7 7 



1- 



Cio 



u/ 2AT t 
lN3C 8 - 
u) 2N C . 

uJ \N C 



Ab 



(/) 



Ab 



(/) 



3 

2N C 



Ab 



if) 



where u = 1 — u. The Ac^ and A&P' denote terms 
depending on a s generated by matching from Hy/, and 
will be considered at 0{a s {mb)) below. The displayed 



( ft 

terms in c\ and b ( ^ ' correspond to (73 4 9,10 terms in 



df) 



the first line of Eq. (|1 1|), and Ac^ and Afey ; include 
the a s C\ 2 80 terms. 



IV. ENDPOINT BEHAVIOR OF (j M (z) 



In order to model the generalized form factors like 
Cj M (z) it is useful to know their behavior in the end- 
points z — > and z — ► 1. This behavior along 
with that of <^> M (it) determines whether the convolutions 

/ duc[ f \u)<t> M {u) and / dudz b[ f) Ju, z)(f M (z)(f> M {u) in 



the factorization theorem in Eq. (42) converge naively or 
require zero-bin subtractions |30|. At tree level c^\u) ~ 



1 since it is independent of u, while b\ '(u,z) 
independent of z (for the electroweak coefficient b^ re 
place the u^ 1 by u" 1 ). At one loop the scaling behavior 
becomes c\ (u) ^(uu)^ 1 and b\ (it, 2;) ~ (zu) _1 as dis- 
cussed in section [V] below. Known two-loop corrections 
do not modify these one- loop scaling results [Q . 

Using a factorization of the generalized form factor we 
can connect the scaling of Cj M to that of (j> M (it). Sepa- 
rating the scales A 2 <C m^A gives the factorization theo- 



j-BM 



(^) 



d.T / dfc J(z,a;, fc_|_) 
^0 



i B (fc + )0 Ml (^), 
(46) 



where J is a "jet function" that can be determined as a 
power series in a s ([ii) where /ij ~ \/mbh.. At tree level 
J is given by 



J{z, x, k+) = 5(x — z)Ci 



a s (M) 7r 
N c zk+ 



Using this result in Eq. (]Am gives 



47ra s (/it) 



3mR z 



(47) 



(48) 



which demonstrates that that the endpoint scaling 
4> M {z) ~ (zz) implies Cj M ( z ) ~ z - Beyond tree level the 
scaling is determined by corrections to J(z, x,k+), which 
is currently known to one-loop order J40|, [41[ and still 
yields (f M (z) ~ z. This scaling is expected to persist to 
all orders. Evidence for this comes from the argument of 
Ref. p5] that is based on an assumed correspondence be- 
tween soft and collincar endpoint singularities in the form 
factor. A strong argument for why these corresponding 
contributions must always arise when endpoint singulari- 
ties appear was given in Ref. [ pp| . Endpoint singularities 
are simply an artifact of not properly separating momen- 
tum regions in the effective theory, and arise in situations 
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where a collinear momentum generates a double count- 
ing with modes that account for the region where the 
momentum is soft (and vice versa). In the effective field 
theory this is avoided by including zero-bin subtractions. 
Thus zero-bin subtractions are not expected to arise for 

M (z), whereas the analogous factorization theorem for 
C^ BM to Eq. (ff6|), which exhibits endpoint singularities, 
requires zero-bin subtractions p0|| . 

In Eq. (^6|) the normalization depends on decay con- 
stants and the inverse moment parameter 



leaving the details to follow in the next two sections. For 
convenience we will use the following decomposition of 
the a s corrections 



Ab 



Ac 
Ab 



(lc) 



Ac { t ] 



(lc) 



Ab 



(lu) 



Ac 
A6 



(is) 

4 

(Iff) 



0{c? s ), 
0{al). 



(50) 



dk + <Mfc + )/(3fc+) = 1/(3A B ) ~ A 



QCD 



(49) 



Here superscripts (lc) and (lu) denote the one loop con- 
tribution due to charm and up quark loops respectively, 
while (lg) refers to 0{a s ) corrections due to operator 
0« . We summarize the results for these terms in the 



A, 



67V c 7r y u 



Unfortunately use of Eq. @ to determine the normal- NDR and HV schemes. For Ac[ f) we have 
ization, Q" , has a large uncertainty due to the unknown 
parameter /3b and the a s (pi) expansion. The pi depen- 
dence can be reduced by using one-loop results for the 
jet function, but this introduces additional uncertainty 
from other moments of 0s(fc + ). 

In this paper we avoid using Eq. (|i^) and instead re- 
duce the hadronic uncertainties by dealing directly with 
the normalization parameters £ BM and Cj M determined 
from data. The analysis in this section constrains the 
model that we construct for the shape of (f M (z), as dis- 
cussed in section XI below. 



Ac[ lu) 



S, 



Ac^ 



6iV c 7r [ * y ' ' 3 
C 8g m b C F a s (p) 



2N r TTU 



(51) 



V. SUMMARY OF LEADING-ORDER 
ONE-LOOP COEFFICIENTS 

In this section we summarize the main results of our com- 
putation of the fourth term on line one of Eq. ( |ll|) , while 



which agrees with the computation in Ref. jl7j and the 
verification in Ref. |l6| . In the NDR scheme the constant 
S^dr = while in the HV scheme Sf v = 1. One of our 
main result is the corresponding corrections for Ab 4 , 



Ab 



(lc) 



A (/) 2i 

c 4tt 



h%(u,z,p)[3C 2 (2z-l)+C 1 (7z+l)} h\{u,z, p) P C 1 (16u+16z-27) h c 2 {u,z,p)C 1 (l-8u) 



9 uz 



27u 2 z 



27 u 



h L 3 (u, z, p){9C 2 uz + C 1 [uz{16z-l) + 2{16z 2 ~25z+27)p}} [27C 2 (2z-l)+C 1 {(6m+55+S b iU+S b2 )z-18}} 



Ab 



(lu) 



54 u 2 (\-z)z 

a s f h%(u,z)C 1 (l-8u) h%(u,z){9C 2 +Ci(16z-l)} 



162 uz 



4tt 



27 u 



Mu{l-z) 



[27C 2 (2z-l) + Ci{(64u+55 + 5 6 ii7+5 62 )z-18}] 



a^ = -(a(/) + aW)) 



162 uz 

Csgfhb C F a s (p) 
mi, 2N c tt u 



r 



(52) 



In the NDR scheme the constants S£ DR = S fc N 2 DR = 



while in the HV scheme Sf^ = -48, and S^ 2 V = 6. The 

contributions from 0$ g , Acjj 1 and Ab^ 9 ^ , are generated 
at tree level and so are scheme independent at this order. 



In Eqs. (H^JH), p = ml/ml and 
ho(u,z,p) 
Li 2 



P_ 

uz 



Li s 



i - g(p/uz\ 



-Li, 



Lio 



1 + g(p/uz)^ 



1 + g{p/u\ 



i-g(p/u) 



hl{u,z,p) = G(uz,p) - In (—2-) 
ho(u, z, p) — G(uz, p) , 



(53) 
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where we have the usual massive and massless fermion 
loop functions 



G(x,p) = \n(^)-26(4p-x)g(p/x)cot- 1 [g(p/x)] 



+ 2-9(x-4p)g(p/x) <j In 
^mix. 



i + g(p/x) 
l-g(p/x) 



2 

G (x) = 2 + lnf^-) + «r, 
\mixJ 



(54) 



with g{x) = yl — Ax and g(x) = \/Ax — 1. The hi func- 
tions are given in terms of loop integrals in Appendix |a|. 
The factors of a s in Eqs. (51 5^) should be evaluated 
at p ~ mf,. One can also look at the endpoint power 
law behavior of these matching coefficients, for which 

Ab^ g) - l/u, and 



(lc),(l«) 



1, Aci ls) 



we find Ac^ 

Ab (lc),(lu) 2 i/^z). 

The SCET Wilson coefficients c{ f) and b[ f) should not 
depend on the 75-scheme choice for Hw- From the point 
of view of the electroweak Hamiltonian, the scheme de- 
pendence in Eqs. (|5l],^2|) corresponds to that in matrix 
elements, and is compensated by scheme dependence of 
the electroweak Hamiltqnian's Wilson coefficients. At 
lowest order, C4 and b^J ' in Eqs. (^4|,|^|) depend on the 

penguin coefficients C3 and C4. Since we are calculating 

( f) if) 

c\ and b\ to order a s we need to take into account 
the scheme-dependence of 6^4 up to order a s , which for 
p ~ nib is given by 



£fNDR 



c 3 HV 



a 



NDR 



= a 



HV 



a s Ci 
36tt 

12tt 



0{a % s ) 



(55) 



Thus at this order the results in Eqs. (fl4|,[l5|) have to be 



used in the same scheme as the Aci^ and Ab^P 



( f) ( f) 

tions. In c\ and b\ we find that it is only the scheme 
independent combinations 



„(/) 



df) 



-a; 



-A 



(/) 



C 4 
C 4 



3 

1+u 
3u 



2a s C 1 S c 



27-k 

a s Ci(S b iu + S b2 ) 



(56) 



648?ru 



that occur. This demonstrates that our final results are 
independent of whether we use the NDR or HV scheme. 



VI. CONTRIBUTIONS OF THE CHARM LOOP 
(NDR SCHEME) 

In this section we present the matching calculation in 
the NDR scheme, where we have an anticommuting 75 in 
d = 4 — 2e dimensions. For Ac\ we calculate the full- 
theory graphs in Fig. [l]a,b,c and match them onto the 
SCET graph in Fig. g. In order to ensure that the NDR 
scheme is consistent, it is important to avoid computing 



a) \ d,s q u 



c) 



d,s 

q 

<&wmm 



b 
q 



1 k q 

■innnr 1 



q 



q 



q 



FIG. 1: Full-theory graphs for the matching onto the short- 
distance penguin coefficient C4 at 0(a s ). a) loop graph with 
charm quarks and or up-quarks and C" 2 , b) counterterm 
graph with Odg, and c) graph with Cs 3 - 



,'q 



d,s I 
b !/ q 



2G F Jd,8) p 
V2 C 4 ^ 



FIG. 2: SCET Feynman rule for the Qj~j operator where the 
spinors are F = [u-^ ^P L u^] {u^ ifiP L v^]. 

traces from a closed fermion loop, tr^^^^^T^l (see 
Ref. |Q for a review). In the basis shown in Eq. (p3|) the 
charm fields are not in the same bilinear, so Fig. Ilia does 
not involve a trace. Since we are treating m c ~ mb, there 
are no corresponding loop graphs in the effective theory. 
Possible loop corrections stemming from NRQCD loops 
vanish at leading order in the power counting, as dis- 
cussed later in section [fx] on long-distance charm contri- 
butions. 

To renormalize Fig. ^a we use the operator 



O 



DG 



J[^,. 9 C7H7,(l-7 5 )&, 



correc- which appears in the electro-weak Hamiltonian as 



G 



H c .t. = -^j= ^2 ^CdgOdg ■ 



(57) 



(58) 



p— u,c 



In the standard basis for Hw given in Eq. (p3| ) the op- 
erator Odg is redundant and has been removed using 
the gluon equation of motion (corresponding to an on- 
shell basis of operators). For our computation we keep 
Odg with a pure counterterm coefficient, SCdg, an d use 
it for rcnormalization. This has the advantage that the 
counterterm graphs maintain their topological correspon- 
dence with the divergent loops (when the divergent loop 
is shrunk to a point). Furthermore, it allows us to obtain 
the desired matching results while avoiding the use of 
rf-dimensional Fierz relations with evanescent operators. 
At the end of the computation we remove Odg following 
Ref. p8| , by writing it in terms of four-quark operators 
and an operator that vanishes by the equations of mo- 
tion, [D^,G^\ = —gT a qj' y T a q, and transforming 
the two-loop anomalous dimension back to that for the 
standard basis. As shown in Ref. Ji3] this gives the usual 
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two- loop anomalous dimension in the NDR scheme [ p9| . 
Thus our NDR scheme coefficients arc the standard ones. 

The graph in Fig. [l]b involves an insertion of the oper- 
ator Odg with counterterm coefficient SCdGi where 



C D a G e O DG ^ ,A°} = C DG 



Orcn 
DG 



sc 



DG 



U DG- 



where definitions for the loop integrals Io(q 2 ) etc. are 
given in the appendix. The explicit 1/e comes from the 
counterterm graph and cancels the divergence in Io(q 2 ). 
In terms of momentum fractions Eq. (|5^) yields the NDR 
result given above in Eq. (|5l|). 



Thus SCdg corresponds to a combination of a countert- 
erm for composite operator renormalization, and wave- 
function renormalization, which for our purposes are not 
required separately. The choice SCdg = — 4Ci/(3(47r) 2 e) 
cancels the 1 /e divergence in Fig. [l]a. The same value for 
SCdg will be used to renormalize the Hw graphs needed 
for the matching computation for Ab± below. 

At this order in a s we only have the tree level graph 
shown in Fig. | on the SCETi side. Matching Fig. | 
and Fig. || gives the following contribution to the SCET 
in the notation of Eq. (Bffl) : 



coefficient 




Aci lc) = - 








Ac^=- 





CfQ-sCi (2m 



6ttN c 



r 



:[I (q 2 )-I ]+I ( q 2 )- — - 



6irN c 



m \9t:u J 



(59) 



( f) 

Next consider the computation of Ab 4 which comes 
from matching the full theory loop graphs in Fig. || and 
counterterm diagrams with Odg shown in Fig. ^, onto 
the tree level SCET graph in Fig. [|. We refer to the 
graphs in Fig. (H) as G a , Gb etc. and those in Fig. (||) as 
SG ab , SG C etc. 



The results for the graphs with quark-loops, Fig. fm-f, are 



G a 
Gi, 

SG ab 

G c + SG C 
Gd + SGd 

G e + SG e 

Gf + SG f 



Gg + SGg 



ig 1 r _ , , 

ig' c 



q 2 (Air) 
7 3 i 
q~ (4tt) 2 
2iCig 3 f abc 
3(47r)Ve 
^2g 3 (2C 1 )f abc 



j. ) ( CiT a T b 
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P L b v 
1 PlK 



3q 2 (p + q) 2 
X {q 2 P a + P -q(p a -q a )} 

ig 3 (2C 1 ) 



[q n liT b q n ] [dn{i - i>)P L T c b v ] 
[dnJ a P L T c b v ][q n jiT b qn] 



2m 2 



6m b u 

«fl 3 (2Ci) 
6mb 

ig 3 (2C 1 ) 

6mb(l— u) 
i9 3 (2C 1 ) 
6m&(l — u) 



[dn4-fi^T a T b P L b v ] [q nl i_„,T b qn 



(An) 2 \(p + q) 2 
2ml 



(I ((p+q) 2 )-Io)+I ((p+q) 2 ) 



1 4 

3 



(An) 2 { q 2 
i \ 2m\ 



(I (q 2 ) - I ) + W) - - - i 



1 4 
3 



[dni^lT^PLW] [q n ^T b qil ] — s 



(7 (g 2 )-/o) + /o(g 2 )-i-^ 



(60) 



[d n iT b P L b v ] [q nl lT b T a q n ] 
[dni,T b P L b v ] [q nl lT a T b q n ] 



2ml 



i)2 -(Io((p + q) 2 ) Io) + h((p+q) 2 ) 



2ml 



(4tt) 2 \ (p + q)- 



;(I ((p + q) 2 ) - I ) + Io((p+q) 2 ) r 
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FIG. 3: Graphs for the matching onto the penguin coefficient 64 at 0(a s ). Not drawn are i) graphs with no gluon attached to 
the quark loop, which vanish in the NDR and HV schemes due to the chirality, and ii) graphs with only the gluon of momentum 
p radiated from the quark loop which also vanish. Here the momentum fraction of the gluon is z with n- p — mtz, the (/-quark 
has fraction 2, with n • q = mtz, the g-quark has momentum fraction u, so n ■ q = rribU, and the d or s quark has fraction it. 
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FIG. 4: Counterterm graphs for Fig. ^ involving the operator Odg- The graph labeled ab) is the counterterm for the sum of 
graphs a) and b) in Fig. 
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In G a and Gb we have objects /" and f% which are defined by 
A^g^c) = |2m c ^ -p __jj 0+ ^_ + _ + _-- §5 _-__j [Jo((p+ 9 ))-/o( 9 )] 



3p-<? 



[/«>(&>+*)*) - Jo] - ^V^o(( P +?) 2 ) + (y(<f -*») - |£) , 



h ^ mc) = - m < (77 J + V + ~W " W " wJ [ o{iP q) ] " 0(9 )] 



2„Q 



2mV 
3p-g 



/o((p+<z) 2 ) - Jo] - ^/^^((P+g) 2 ) + Q(g- - P ") - } • (61) 



Results for the loop integrals Jo and Jo are given in ap- 
pendix |A[ Note that the only contributions from the 8G- 
countcrtcrm graphs are explicit 1/e's, which exactly can- 
cel divergences due to the loop integrals. We have made 
some simplifications in the expressions for the last four 
graphs, where u is the momentum fraction of the d-quark 
in Mi and z is the momentum fraction of the quark in 
M2 ■ The result G g + SG g corresponds to the contribution 
of the expansion of the external full QCD q-quark field 
onto the n-collinear quark field in SCET. For AS = 1 
decays with f = s one makes the replacement dn — > s ri 
in Eq. ©. 

Next we sum the graphs in Eq. ( p0[ ) and Fierz them 
to match onto the SCET operator , and thus obtain 



Ab{ f) . The results for G a + G b + 5G a b, G C + 5G C , etc. are 



and then read off the prcfactor of the spinors P Xa shown 
in Fig. [|to obtain Ab/\ A few useful Fierz formulas are 
listed in appendix |A|. The definitions for all momentum 
fractions are summarized in the caption of Fig. ^. The 
result in the NDR scheme is 



d,s 



f.->-5 = Z^GfsI b (d,s) p 



Xa 



in b 



>1-Z 



FIG. 5: SCET Feynman rule for the Qy] operator where the 



finite as e — > 0, so we can Fierz them in 4-dimensions, spinor product is P* 



[ui d) T^ ± P L ui b) ][ui q) 4PLvi q) ] 



. (ic) = A(/) a. J m 2 c J [3C 2 (2z-l) + C 1 (7z + l)} _ (I (g 2 ) - I Q )pC 1 (16u+16z-27) 
4 c 4tt I 9uz 27 u 2 z 



(W 2 ) -i)Ci(l-8u) , (I ((p + g) 2 ) - Io(q 2 )){9C 2 uz + C 1 [uz(16z-l ) + 2(Wz 2 -25z+27)p}} 

(62) 



27m 54u 2 (1-z)z 
[27C 2 (22-1)+Ci{(64u+55)z-18}]1 



162 uz 

The result for the up-quark penguin loops is simply obtained by taking m c — > which gives 

AI.C1-) - xW as / (J^(g 2 )-7)Ci(l-8") , (4" ) ((P + g) 2 )-4" ) (g 2 )){9g 2 + Ci(16z-l)} 
A&4 ~ K 27S + teu(l-z) 

[27C 2 (2z-l) + Ci{(64u+55)z-18}]l 
H 162 uz J 

These results for Ab£ c) and Ab[ lu) in NDR are summarized in a more compact notation in Eq. (p2|). 



(63) 
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Next we quote the results for contributions of the graphs with the operator Os g , Fig. |||h-m give 



Gh = ffT 6 (4-) r b ndnPRT c b v ][q n liT b qn ] (64) 
%-K l mi \uzJ 



ig 3 C Sg m b ( 1 \ rabc [J ( A u _J>J 1P „ TO u- ur r b r 



l6Tr 2 ml 



II z 



Gi = ~ 9 * C 2 * a ? b (~) [d ri P R T b b v ][qnliT b T« qfl ] 
Gm = ~( C 2 8 f b (~) [d Tl P R T b b v ][q nl lT«T b q Tl ] 
Go = ^^^f 1 --^ [dnlil^P R T a T b b v ][q nlUl T b qfl 



8ir 2 ml \z(z + uz) 
Gk = G n = 



Here mt, is the MS mass which always accompanies Cg g , 
while mb is a short-distance threshold mass (for which 
we will use the lS-mass). For / = s we take d fi — + Sr. 
Fierzing the results in Eq. ( |64] ) and matching onto Fig. || 
we obtain the contribution to Ab^ due to 0% g , 

A^ = -(A(')+A(/>)^(f^V (65) 

TOf, \ 97TU / 

which turns out to be identical to the matching result for 
Aci l9) . 

VII. CONTRIBUTIONS OF THE CHARM LOOP 
(HV SCHEME) 

Next we repeat the calculation of the previous section 
using the HV scheme. In the HV scheme, 75 anticom- 
mutes with Dirac matrices in 4-dimcnsions, and com- 
mutes with the Dirac matrices in the remaining (— 2e)- 
dimensions. Here we can consistently handle traces with 
7 5 in d 7^ 4 dimensions. In the HV scheme the Dirac 
matrices in the (V — A) interactions in the weak Hamil- 
tonian are taken in four-dimensions, while all 7-matrices 
from the QCD and QED Lagrangians are in ci-dimensions 



(see Ref. ]47| for a review of the HV scheme). We will 
perform the computation in two operator bases, namely 
the original one 0\ 2 in Eq. (^3|) and a different basis 
given by 

Of = \db\v-A\w\v-A, 

61 = [dT a b] v _ A \pT a p] v _ A , (66) 

with Wilson coefficients G\ and C2. Notice that in 
addition to the different Fierz order, the <D\ 2 basis is 
also color Fierzed with respect to the 0\ 2 basis. In 

the 0\ 2 basis we need to evaluate fermion traces like 
tr[7 AI 7 l/ 7 a 7 /3 7 5 ]. As we explain below, the answer for 
the matching computation in either of these bases is the 
same in the HV-scheme. Renormalization with the op- 
erator Odg goes the same way as in the NDR scheme 
except that we replace C\ = C2/2. 

The result for the charm and up-quark loop in the HV 
scheme is well known, see e.g. Rcf. J39|, and with either 
basis gives the same result for Ac^ 1 ^ and Ac^ 1 "-* as the 
NDR scheme but with Sf v = +1 in Eq. (|l|). 

For A64 we start with the computation in the 0\ 2 
basis. The graphs in Figs. |3| and |] give 
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G a + Gb + SG a b 
2iC x g 3 a 



Xp,p 



(47T) 2 



(Pp - Qp) 



(ml Jo + -) - (q p +p p ) 7 ^ z (I ((p + q) 2 ) - Io{q 2 )) 
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[dnJaPLbv] [qnl±pT a q n ] 



^29 ), a bc c a\np 



(4tt) 2 i 
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^-^-(mt Jo+l)~ (qp+Pp)^-(Io((p + q) 2 ) Io(q 2 )) 



;(p + ? J + q - q - 
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2m 2 c q 2 



[dnj a P L T c b v ] [q n J± fJ ,T b q n 



4p.q 



'ip.q 6p.q 3(p.q) 2 \2(p.q) 2 



p a Ul ((p + qy)-Io(q 2 )) 



1 \ 2m 2 c p° 



op.q 18 



ip.q 

[dnJ a P L T C b v ] [qnllT b qfl 



(h((p+q) 2 ) - lo) - ( m 2 c (p a - q a ) + ^P° ) J> 



p.q 



(67) 



where we use e 0123 = +1. In Figs, ^c-g only 0\ contributes, and we have 
-2C 2 g 3 f abc 



G c + SG C 

G d + 5G d 
G e + SG e 
G f + SG f 



Gg + SGg 



3q 2 (p + q) 2 
X {q 2 p a +p-q( P a -q a )} 
*C 2 g 3 

6m,bU 

iC2 9 3 
6mb 



dn"f a PLT c b v ] [q n J±T b q fl ] 



(4^{(^ (/0(b+<Z)2K/ ° ) + /0((P+,Z)2K ^3f (,,S1 



[dni^lT b T a P L b v \ [q n l±pT b q n ] 



^{^(W)-7 ) + / (^)-i-i 



iC f&u) [^T b P L b v ] [q nl iT b T«q,] ^ {j^(Io(( P + q) 2 ) h) + I ((p+q) 2 ) - - 4i 



6m&(l — u) 
iC 2 g 



6mb(l — u) 



^— [d fl 4T b P L b v ] [q nl lT a T b q n ] — 

1— Ml L J L J 47T) 



2m 



(4tt) 2 \ {p + q)- 



e 3 

;(I ((p + q) 2 ) - Io) + Io(( P +q) 2 ) - - - \ 

e 3 



r 



Results for Jo and Jo can be found in appendix |X| Again 
we are free to Fierz these finite results in 4-dimensions. 
Computing Ac± ' and Ab^ from these expressions gives 
the results summarized in Eq. (p2|). 

Alternatively one can do the HV scheme calculation 
in the same basis 0\ 2 as the NDR scheme calcula- 
tion. Although there are no fermion loops in this basis 
the HV scheme computation does differ from the NDR 
scheme. For each graph the results differ due to an extra 
0(e) term generated in manipulating the Dirac matrices. 
Therefore it is easy to quote the HV scheme results ob- 
tained in this basis, as replacements to be made in the 
in NDR result. For /f in Eq. (|6|) we should replace 



For graphs G c + SG C to G g + 5G g in Eq. ( |60| ) and for 

A, 



in Eq. J59| ) we replace 
4 

~ 3 



(71) 



Finally the HV scheme result for Ab^ p ^ in Eqns. (|6 
and (K33) is obtained by the replacement 



[27C 2 (2z-l) + Ci{(64u+55)z-18}] 
162 uz 

[27C 2 (2z-l) + Ci{(16£t+61)z-18}] 
162 



(72) 



This is same as the result that we obtained from the HV 
scheme calculation in the 0\ 2 basis. As discussed earlier 

in section yA the scheme dependence in C3 and C4 which 

( f ) 

(69) appear at tree level accounts for the shifts in Ac 4 and 
and for the HV scheme result is obtained by replacing Ab{ f) given by Eqs. @ and @, thus leaving the SCET 



(q a ~p a ) 



q 2 p a 
6p-q 



(q a -p a ) 



q 2 p a 
6p-q 



q 2 p a 
6p-q 



(70) 



Wilson coefficients independent of choice of NDR or HV 
scheme. The calculation in the HV scheme in the 0\ 2 
basis differs from that in the NDR scheme, and provides 
a non-trivial cross-check on our results. 
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VIII. CHIRALY ENHANCED PENGUINS 

It is well known that certain power corrections have 
the potential to be numerically enhanced in penguin am- 
plitudes. In particular the so-called chiraly enhanced 
terms [ )l7j] , which are formally down by a factor of A/ rrib, 
but are numerically of order /ip/m;,. 2 For the pion 
/i,r(2GeV) = 1.7 GcV, and this can be understood from 
the fact that oc A x rather than Aqcd ; where A x is the 
scale of chiral symmetry breaking. Thus relative to the 
other power corrections these terms have the possibility 
of being magnified by a numerical factor of n^/A ~ 3 — 4. 
A valid factorization theorem for the complete set of chi- 
raly enhanced corrections has not yet been derived, be- 
cause previous attempts encountered endpoint singulari- 
ties |l7j . In this section we derive a factorization theorem 
for chiraly enhanced tree and penguin amplitudes that 
does not suffer from endpoint singularities. Our analy- 
sis uses factorization in SCETt and the complete result 
involves only one additional generalized form factor and 
one light-cone meson distribution beyond those occurring 
at leading order. 

To consider chiraly enhanced operators in SCET we 
can work with a complete basis of operators suppressed 
by one power of A/m b , and then look for all operators 
with a If i_ in the light-quark bilinear as explained in 
Ref. p9| . This provides a unique way to determine the 
contributions that are chiraly enhanced, without invok- 
ing the Wandzura-Wilczek (WW) approximation as 
was done in Ref. 

We therefore construct a complete basis of operators 
with one P", starting with the field structures: 



(73) 



Only the color structures shown are required at this or- 
der. Operators with a igB^j_ are needed at the same 
order as operators without, because of the additional 
suppression of the non-2?"j^ operators in the matrix el- 
ement of the required time-ordered products. This is 
the same situation which we described already at lead- 
ing order in A/m b in Eq. (||). Note that in Eq. (|7|) 
we do not consider other operators with V± or d± since 
they are not chirally enhanced. To perform the match- 
ing we work with a basis of four-quark operators of def- 
inite chirality, where the possibilities are inherited from 
the full electroweak Hamiltonian: (LH)(LL), (LH)(RR), 
or (RH)(LR). Here the order corresponds to the quark 
fields in Eq. (J73j) and we do not assign a chirality to the 



2 Although chiraly enhanced penguin contributions are large, for 
tree amplitudes they are numerically the same size as other ex- 
pected power corrections as emphasized in Ref. 



heavy quark denoted by H. With definite chirality a 
complete basis of Dirac structures includes 



r„e{l,7^}, I\e{^7±}, 



(74) 



where r ft = r/i contributes only to (LH)(LL) and 
(LH)(RR), while T n = rij 1 ^ contributes only to 
(RH)(LR). 

First let's construct a complete basis of the Q^^-type 
operators in Eq. (|73]). Here T n (g) T fi must have a _L (i 
index, and we find the basis 









- O ilx) -e 


o {lx) 

^3(qfq) 






rrib 



'_L Snt^ ] ' 



(75) 



where we have f = d,s. The {qfq) subscripts on the op- 
erators indicate the flavors of the light quarks, and the 



basis has in addition the operators Q 



3(/L) and Qafrfvy 
Whenever a flavor label q appears we implicitly sum over 

q = u,d,s. The operators give contributions to 

PP, PV , and VqVq final states, whereas Q y 3 only con- 
tribute for transversely polarized vector mesons. If oper- 
ators that produce n-isosinglet mesons are included, we 



have in addition Q^ X f qq ) and Qftf qq y 



Since the Q^a 



operators have right-handed quarks, only 05,6.7,8 m the 
electroweak Hamiltonian can contribute to their match- 
ing at tree level, while other operators start contributing 
at one-loop. 

Next we construct a complete basis for the Q^i* 1 an d 
Qg2 ''"type operators. For chirality (LH)(LL) we must 
have r„ = 1 and Tn = 7/ and we have two choices for 
contracting the _L indices, or ie°f . To avoid having 

the epsilon symbol in our basis we trade for a pair 
of 7j_'s. Here the possible flavor structures are (qfq), 
(ufu), (fuu), and (fqq) from matching the operators in 
the original Hyy- For chirality (LH)(RR) the same Dirac 
basis applies, with flavor choices (fuu) and (fqq). The 
latter flavor structure only produces n-collinear isosinglct 
mesons. Finally for (RH)(LR) we must have r„ = 7^ 
and T n = 7/7^ and there arc only two inequivalcnt ways 
of contracting the _L indices (a(3uv) . This follows since 
contractions with an ie± do not lead to independent 
structures because of the fixed chirality, and the identity 
7^7^ Pl <g> iftrf-hPn = which allows an additional con- 
traction to be eliminated. For (RH)(LR) only the flavor 
structure (qfq) contributes. All together these results 
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lead us to define the basis 
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(76) 



plus operators with the same Dirac structure but differ- 
ent flavors, Qi(^^and Q^ u ) for i = 1 and 5. The 
operators in Eq. (|7q) also incorporate electroweak pen- 
guins, since we can write e q qq = uu — 1/3 qq. Oper- 
ators Q [2x) A , contribute for B -> PP, B -> PV, and 

— » VoVb decays, whereas the operators only 
contribute for decays with transverse vectors in the fi- 
nal state, B — > VtVt- If n-isosinglet operators are in- 
cluded we have in addition the operators Q^j^ qq ) where 
% = 1,2,5,6. 

For the basis in Eq. ( |76| ) we have only written operators 
that contribute to B decays. The remaining operators 
which only contribute for weak B* decays are 



^Hqfq) 
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L 1 
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2(fuu) m 
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+ i [^^i ^^^] [AL^SsU] . (77) 



to non-isosinglet final states. We demonstrate that the 

Q\(p{ only contribute for B* decays in appendix |B|. 

The Hamiltonian for the full basis of (1%) and (2x) 
type-operators contributing to B-dccays is 



V2 



[dWl,2,3] c f(F)K')^(F)K ) 



(78) 



„x 

-i(F) 



J i(F) 



where the indices run over the operator number i and 
possibilities for the flavors F for the Q^f)' s shown in 
Eqs. @ and ©, and , " 
Wilson coefficients. 

Next we match from Hw onto the operators in 
Eqs. ( [75| ) and ( |76] ) to determine the Wilson coefficients 
cffp\ and bf,ps. at lowest order in the a s (mb) expansion. 
At lowest order c^. ^ ^ are simply given by the matrix el- 



i(F) 

ements of the Oi 's expanded to next-to-leading order in 
the A power counting with 



?= 1 



1 



(79) 



For the (LH)(LL) and (LH)(RR) chirality only the ex- 
pansion of the fi-bilinear contributes, and for the non- 
isosinglet operators we find 
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(80) 



^i(F) terms denote perturbative correc- 
tions. Numerically they will not always be suppressed 
due to the competition between 63,4 and a s (mb)Ci. For 

the operators Q^^F) on ^ ^5-8 from Hw contribute at 
tree-level since the operator involves right handed quarks. 
We find 
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(81) 



and Q^slqfq) — l e q Q 3^7(9/9)- Taken together the results 
in Eqs. (^) and fl77| ) form a complete basis for decays 



Wc find that the loop and magnetic penguin graphs in 
Fig. [TJ can only contribute to the matching when a factor 
of _L-momcntum is generated by expanding the spinors. 
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FIG. 6: Full theory graphs for the matching onto the short- 
distance penguin coefficients for the Qbi and Qb2 type oper- 
ators. 



They give the following 0(a s ) corrections to the match- 
ing 



b), and d). We find 
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(83) 



and 
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(84) 



where the other Ac^ F ) coefficients are zero at this level, 
and S c = for the NDR scheme while 5 C = 1 for the HV 
scheme. The scheme dependence in these results cancels 
against that in the tree level 63,4 terms in <^i q fq\ just as 

for the LO c\ Wilson coefficient, and also in an identical 
manner against the scheme dependence in the tree-level 



C5 6 terms in c 



Results for the Wilson coefficients for cases with isosin- 
glet n-mesons can be determined in an analogous way, 
but in this case operators with gluons also become nec- 
essary. We leave this for future work. 

Since we are only using factorization of effects at 
u nib the matrix elements of the operators Q^pl and 



Qi(F) §* ve hadronic parameters. For the 77,-collinear field 
we need the matrix elements in Eq. (flO]), and the cor- 
responding results with right-handed light quarks and 
Clebsch-Gordan coefficients 



1(9/9)- 



Next consider the matching calculation which deter- 
mines b x ^ F y At tree-level this involves computing the 
graphs in Fig. ^, and involves non-zero contributions from 
expanding the propagators in graphs a), b), and c), and 
from expanding the spinors with Eq. (|79|) for graphs a), 



(M n \Ti [q^pv] \B) = C*f5 Ul m B ( BM , (85) 



We also need a new form factor ( BM (z) and the chiral- 
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TABLE I: Hard functions for the chiraly enhanced amplitudes in Eq. (^) for B° and B decays to PP and PI 7 channels. We 
have not listed results for R{ 2 , but they have the same Clebsch-Gordan coefficients as i?i,2 and so can be simply obtained by 



the replacements c^iqfq) 



& 3(9/9) aIld C 2(9/9) 



^4(9/9) * n * ne commns above. 



enhanced function <^£, (x) defined by 



PV channels 
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(P(P)\ QL 2 ^iC 1°) = —* C lm' 5*** /pMp < P (^), 



where <5 WW ' = <5(cij— cj'— mg) and the momentum fractions 
z = uJi/n-pM and x = 072/?^ ■ VP- where C^ L M and C^ q , 
are the same Clebsch-Gordan coefficients that appeared 
already at leading order in Eq. (|4^). For the chiral- 
enhanced distribution function 0^ we used the definition 
in Ref. (29) . and take the other twist-3 meson distribution 
to be the three-body (j> 3M (x,x) which does not generate 
chiraly enhanced contributions, fn a more traditional 
basis there is a redundancy at this order in f/mf, (see 
for example Q), and 4>p P (x) = 3x[0p (x) + cj)^ / (x)/6 + 

2/3p/(/pMp) / dy' '/y'fap{y - y',y)]- In the Wandzura- 
Wilczck approximation one would set (f> pp (x) = 6x(i — x). 

Taking the matrix element of the above operators leads 
to a factorization theorem for the chiraly enhanced am- 
plitude for non-isosinglet charmless -B-decays to PP and 



A X (B -> MxM 2 ) = 



'Fm B 
V2 



(87) 



JduR 1 (u)^{u) + {l^2) 



dudz i?i 7 (w,z)Cj M2 (z)^ 1 (u) + (l 2) 



Stub 

BMllML f du dz R x (u, z)C^f M2 {z)4> Ml («) + (1^2) 
6m B J x 



This amplitude only includes the chiraly-enhanccd power 
corrections where factors of \im are generated by pscudo- 
scalars, and so for vectors we define piy = 0. (Note that 
we include the symmetry factor of 1/2 in the branching 
ratio prefactor for B — > 7r°7r° rather than the amplitude.) 
In terms of Clebsch-Gordan coefficients for the different 
final states, the hard functions Ri, i?/,and Rj for the 
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chiraly enhanced amplitudes are 
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(88) 



Summation over g = u, d, s is implicit. Results for these 
hard functions in different channels are listed in Table |. 
Equation (^7|) with Eq. (^) corresponds to the contri- 
butions given in the second line of Eq. (O) (when we 

( ft — 
extract the coefficients of the Ac terms). 

From the matching results we find that the endpoint 
behavior of the Wilson coefficients is ~ 1/ (uu) and 
bf/ F \ ~ \/{zuu). Since we know the endpoint behav- 
ior (j) pp {u) ~ uu and 0( z ) ~ 2, it remains to determine 
the behavior of ( x (z). The operator defining C x (z) has 
an extra V±/n ■ V relative to the operator defining the 
distribution Qj(z). Now from the collinear power count- 
ing V± <C n ■ V, so consistency of the power counting 
in SCETi implies that the scaling of ( x (z) as z — > 
and z -» 1 can be no worse than C,j(z). Thus we take 
Cx( z ) ~ z - This demonstrates that all the terms in the 
factorization theorem for chiraly enhanced penguin and 
tree contributions given in Eq. ( |87| ) converge, just like 
the leading order factorization theorem in Eq. (p2|). In 
appendix we argue that the same conclusion about the 
z-convolution is obtained if one considers the direct com- 
putation of Cx(z) m SCETn. 

As already noted, the operators in Eqs. (75) and (|7q ) 
also generate contributions with two transverse vectors 
in the final state. To take the matrix element of these 
terms requires 



(89) 



where u = and three form factors 



(V\qL^lbv\B) = C^S^mB Cf V 'ef ± , (90) 



{V\q^igBi 



b v \B) = -CS UlUt m B Cj±(*)e±. 



where z = Ui/mb- Thus, our complete basis of operators 
with V a terms generates a contribution to the amplitude 
to produce two transverse vector mesons that involve two 
types of light-cone meson distributions, and three types 
of form-factors. These analogs of the chiraly enhanced 
terms were displayed as the contributions on the second 
line of Eq. (jl2|). Our analysis demonstrates that only 
these terms will be generated from the V± operators con- 
sidered in this section, however a full analysis of these 
terms will not be given here. Hence we have not both- 
ered to specify the Clebsch-Gordan coefficients C and C 
relative to our other conventions. 



IX. LONG-DISTANCE CHARM 

In order to properly determine the short-distance coef- 
ficients by matching we must make sure that we subtract 
any effective theory diagrams. Earlier we stated that 
there were no SCET loop graphs to subtract. In this 
section we further justify this claim and discuss long- 
distance charm contributions. We take m c ~ nit, and 
so do not have collinear charm quarks. Furthermore, 
graphs with collinear or soft up quarks are power sup- 
pressed. The only remaining term to consider are soft 
non-relativistic charm that propagate in the EFT. While 
a factorization theorem for this type of long distance 
charm effect has not yet been derived, we may never- 
theless match systematically by including in the effective 
theory the proper SCET-NRQCD hybrid operators as 
discussed in |^4|. 

We begin by showing that non-zero contributions from 
the hybrid operators requires a non-zero residual momen- 
tum and hence do not affect the matching computations 
in earlier sections. We may write the momenta of the 
charm quarks as 



u 

Pi 



~M — 1 TV- J. 

P2 — ^ * ' 



(91) 



where q M is the total momentum of the charm quark pair 
and r l is the relative 3-momcntum in the cc rest frame. 
L v ' (q) is the Lorentz boost from the center-of-mass frame 
to the B rest frame and has components 



2m r 



3l 

2m r 



(92) 



When matching onto NRQCD at lowest order in a s we 
generate a generic set of the four quark operators. As 
with the operators in Eqs. ( |35|j3^ ) there will be a set of 
operators with and without gluon external lines. Thus 
there will be two generic forms of the operators 



<l a 

w prod 

n b 

w prod 



(V L{T NR ) x )(In^B^T hl b v ) , (93) 



rVflj T/iz are the possible bilinear Dirac structures for 
NRQCD and a heavy-light bilinear in SCET. The gluon 
field B± has a four vector index that is either contracted 
with Thi or L(Tn R ). For O a the only possible structure 
is 



L(T NR )®T hl = a l Ll<E>^P L 

For the O b operators the possible structures are 
L(T NR ) <g> r w = {1 ® ^P L ,a l Ll <g> faP L } 



(94) 



(95) 
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In addition we have four quark operators that arc gener- 
ated by integrating out one hard gluon exchange. They 
have the general form 



O 



a 4m^ ^ t£ ( rjVR ) X )( g "'"i r "" g "^3) : 



(96) 



For the 0% nn operators the possible structures are 

L(T NR ) <g> r nfi = {1 ® 1, a*L* ® 1, ^4 (» (97) 

while for the 0„ nn operators the possible structures are 

L(T NR )®T nfi = {l^yL'gi^yi^S^}. (98) 

In general both 1 ® 1 and T ® T color structures are 
allowed in the operators O^ od and 0"^ b n . 

A cE then follows from the time ordered products of the 
form 

Tf = 

f^zd*yd*y'T[O^ od (0)O^ n (z) i£M q (y) iC^Jj/)] 

+ Jd 4 zd 4 yd 4 y'T[O^ od (0)O a ann ^) ^ q (v) i^cjiv')}} 
+ Jd 4 zd 4 yT[O a prod (0)O a ann (z),iC ( ^\y)], 

Tf = Jd 4 zd%T[Ol ro MO a ann (z),iC^ q (y)] 

+ T[O^ od (0)O b ann (z),iCf] q (y)]. (99) 

These operators could be factored into soft and collincar 
components, however the details of this factorization will 
not be carried out here. The factorization for semi- 
inclusive decays was discussed in Rcf . . 

Now let us review some aspects of the power counting 
for these terms (refering to the appendix of Ref. [g4j for 
further details). First note that implicit in these opera- 
tors is a label conserving delta function. Recall that the 
NRQCD fields have two large labels |32| Xm,mv which 
have been suppressed in these operators. For instance, 
repristinating the momentum conserving delta function 
and momentum labels for O prodl we write 



prod 



8{m h —loi- n M L^(2m c )). (100) 



Note that these delta functions do not imply that we are 
only including a single point in the phase space, since the 
residual momenta of the HQET and SCET field in the 
operator may flow through the charm loop. Furthermore 
since the residual momentum scales as A ~ tow 2 , the 
fluctuations in the external momenta effectively smear 
over the non-rclativistic region. 

The delta function constraint simplifies the matching 
since the contribution of these hybrid operators to the 




FIG. 7: Charm loop contribution from the non-relativistic 
region. 

matching vanishes at the lowest order in v. To see this, 
we may work slightly away from threshold, by giving the 
heavy quark (without loss of generality) some small resid- 
ual momentum k ~ A , such that the invariant mass of 
the charm quark pair is 



q 2 = (rribV + k — p) 2 = 4to 2 — 



4m| 



i-k . 



(101) 



where n ■ k ~ A and p is the momentum of the d, s quark. 
In deriving Eq. (101) we used the heavy quark equation of 
motion, v-k = 0, and expanded about the threshold, to^— 
n ■ p ~ Am 2 / nib- An explicit calculation of the one loop 
diagram shown in Fig. ^, shows that this contribution is 
proportional to 



\Jq 2 - 4m 2 ~ 



(102) 



Given our scaling, mv 2 ~ A we find that this contribution 
is order v as anticipated by power counting arguments. 
Hence, only if we were interested in matching explicitly 
onto w-suppressed corrections would be need to include 
these hybrid operator diagrams. Since the suppression 
by v leaves these terms larger than other A/mb power 
corrections it is quite feasible that long-distance charm 
contributions are numerically relevant. 

Finally, we explain why in general we expect these long 
distance charm contributions to be complex. As was 
shown in Ref. |38| a sufficient condition for the gener- 
ation of a complex phase is the presence of soft Wilson 
lines in both the n and n directions. In most obscrv- 
ables these Wilson lines cancel, however, it is clear that 
this will not be the case for the long-distance charm con- 
tribution. The underlying reason for the lack of such a 
cancellation is the fact that the charm quark propagate 
over long distances. Thus when we rescale the light quark 
fields 

q n (x) -> Y n {x)q n (x) , q n (y) -> Y n {y)q n (y) , (103) 

the argument of the Wilson lines will be at different po- 
sitions. Furthermore in general the charmed quarks will 
not decouple from the B meson, thus the soft matrix 
element will be of the form 

(0\[vh}(0)[x^}(y)[qYn}(z-)[YX}(0) 
Y n (y)Y 7 t (z-)Y n (0)rJ (y)\B), (104) 

where the spin contractions and the color contractions of 
the fields and Wilson lines have been suppressed. Note 
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that since the charm quark pair propagates over long 
distances the point y^ is displaced away from zero along 
the light cone as well as transverse directions. Since this 
matrix element for long-distance charm knows about the 
two final state hadrons, it can have a nonpcrturbative 
complex phase. 



X. PENGUIN ANNIHILATION 

In this section we review the penguin annihilation con- 
tributions occurring on the third lines of Eqs. ( pd| ) and 

( |l2| ) . For our purposes they are defined as the Ac terms 
in the amplitude obtained when the spectator quark is 
annihilated by the weak operator. These contributions 
start at 0(A/mb) relative to the leading power penguin 
terms. They include the well known terms fB<fi M 4' M 
from spectator annihilation with a subsequent pair cre- 
ation |^5|, |27j], as well as terms of the same parametric 
size /b</>^</> 3M m where the spectator emits an ener- 

The 



getic collinear gluon prior to its annihilation |28 



former require zerobin subtractions to obtain finite con- 
volutions p9t , while the latter do not. With these sub- 
tractions the leading penguin annihilation contribution 
to the amplitude is real. The scheme dependence of 
the zerobin procedure is compensated by terms involv- 
ing the exchange of a soft quark in the annihilation pro- 
cess, Aj^ , which come from time-ordered products in 
SCETi. These time-ordered product terms start at order 
a 2 / nib and have a nonpcrturbative strong rescattering 
phase. At (D(A 2 /m 2 ) one also has chiral-enhanced pen- 



guin annihilation terms fBH-M4>pp4> M 
be factorized with zerobin subtractions as in Ref. 

lM ±m' 



7 which can also 



The annihilation terms /sf f and fBHM<r pl 
do not involve a hard-collincar propagator and so appear 
to be insensitive to the intermediate scale ~ \/nibA. 
However the zero-bin subtraction procedure is needed to 
distinguish soft and collinear regions, and so they are 
not defined independent of A^l nn at 0{a 2 {fXi)). Since 
this procedure has an a s (u,i) expansion we consider all 
penguin annihilation contributions with an expansion at 
the intermediate scale, unlike our analysis of amplitudes 
in earlier sections. This increases the theoretical uncer- 
tainty, and will be accounted for in our error analysis. 

From Rcfs. [^9| and |^8| the penguin annihilation am- 
plitudes at 0{a s ) for B -> M 1 M 2 are 



A 



(i) _ 



/b/Mi/m 2 



(105) 



hcann 



hM 1 fM 2 {H^} 1 M2 ) 
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Here we have the inverse moment of the B light-cone 
distribution 



Pb = 



'dk 



(106) 



M t M 2 



and the five factors (H^ M ), (H^ 12 ), (H hcl 
(ff^ lM2 ), and (H^ lAh ) arc linear combinations of the 
moment parameters 

aM\M; 
r J ic 



Vhcl,3 



oM-l M 2 
Phc2A 



oMi M 2 



f dxdy (a lc + a l+4 ){x, y) 4> Ml (y)(j) M2 (x) , 
Jo 

d x dydy ah f™U 3M Hy,y)<f> Ah (x) 



dxdxdy 
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1—x — x 



3M 2 , 



-^dxdyalr J {x,y)^{y)cj ) M2 {x). 



ft 



Mi M 2 
X2,6 



6 



dxdyal 6 (x,y)<f> Ah (y)cj ) ^(x), (107) 



respectively. The appropriate linear combinations for 
each channel are given by the entries in Tables II, III, 
IV, and V of Ref. 



for A 

(i) 

hcann 



(i) 



Lann and A LaL ™d in Ta- 

. The distribution functions 



ble I of Ref. |2J for A 

(f> M and <j£ * appearing in Eq. (107) were defined above 
in Eqs. ( f40| ) and (|86]), while the three-body distribution 
is defined by the matrix clement 



LJ 3 



hp 4>3p(y,y), 



(108) 



where 5^ U3 = 8(u>i + uj 3 - ui 2 - m b ), y 



uii/rrib, and 

y = —(jjz /rrib. The circle on some of the integrations in 
Eq. (107) indicates the terms which require zero-bin sub- 
tractions. These subtractions modify the hadronic distri- 
butions by inducing dependence on a rapidity parameter, 
which increases the uncertainty from these terms. For our 
numerical analysis we adopt the models used in Refs. |29[] 
and Eq] to determine the /3's in section XII below. 



XI. 



ANALYSIS STRATEGY AND MODELS FOR 
THE SHAPE OF (f M (z) AND C* M («) 



To make predictions for the p M i M ? penguins at lead- 
ing order we need values for the twist- 2 meson distribu- 
tion 4>m(u) and the form factors C BM > Cj M ( z )- To com- 
pute the A/ mb suppressed chiral-enhanced amplitudes we 
need in addition the twist-3 distribution (f> AI p (u) and the 
form factor C£ M (z). 

A common model for (j> M and capturing the essen- 
tial features is given by the first few terms in the Gegen- 
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baucr scries 



using the parameterization 



(j) M {x) = 6x(l-x){l + af (6x-3) + 6af (l-5x+5x 2 )} 
+ 15af (l-14x+56x 2 -84x 3 +42x 4 )}, 

4>f p (x) = 6x(l-x)« p + < p (6x-3) 

+ 6a^ p (l-5x+5x 2 )} , (109) 

where x corresponds to the momentum fraction for the 
quark field (dressed by a Wilson line). Our ranges for 
the model parameters af 1 and af pp are summarized in 
the next section. Note that we include both af 1 and a^f 
in <fi M (x). This allows us to account for data on (x -1 )^ 
which constrains a* + a\ , while varying a * — a\ to obtain 
a range of models. From charge conjugation and isospin, 
^•p(x) = 4>*' p (l - x). Thus we will set a\ = a{ = 0. 
With SU(3) flavor symmetry af = 0, so smaller val- 
ues are adopted for this parameter than are used for af 
(we keep a non-zero af 1 for SU(3) violation, but do not 
include a non-zero a?f or any other odd-moment param- 
eter). These restrictions from charge conjugation also 
apply to the distribution 6^1. To see this we follow the 
same argument given in Ref. |p3| but for the matrix ele- 
ment defining the chiral-enhanccd distribution: 

(TT°\c^c^w n c^c$ l5 f ± s(Lo - r + )w^ n c^c\o) 

= (+l)(7r°|(CW^„) T $75lM(w - P+)(UW n C) T \Q) 
= +(Tr°\tn,p>W n ^ 5 T±6(u + P + )Wte ntP \0) . (110) 



Here V+ = + V . C is the charge-conjugation opera- 
tor and C is the charge-conjugation matrix. To obtain 
the last line we note that the sign from P"£ n = —S.n'P'x* 
cancels the sign from (7„ ) T C = 7:7. Hence the matrix 
element is even under to = (l — 2x)h-p 7r — > —u, which im- 



plies 4>i P (x) = 4>* p (l-x 
and small values of af pp . 

,M 



Hence we have 1 



- n p - 
L lpp ~ 







Note that in Wandzura-Wilczek 
approximation afipp = 1 and af^pp = 0. 

Potentially larger uncertainty comes from the values 
for ( BM and C,f M {z). At lowest order the bi(u,z) coef- 
ficients in Eq. (^2|) are independent of z and thus only 
£ SM and the zeroth z-moment, M , from Eq. (Q) are 
required. This yields two form factor parameters to be 
fit to tree amplitude data. For convolutions with the 
a s Abi(u,z) terms we need more information about the 
z-dependence. However it is still very useful to fit the 
norm Cj M to the nonleptonic tree amplitude data. In 
particular, by only modeling the shape of Cj( z ) we re ~ 
duce the model uncertainty considerably. The fit to non- 
leptonic decay data currently provides the most accurate 
way of determining the normalization of £j M (z). Thus, 
our strategy allows us to make predictions for the pen- 
guin amplitudes based on input about the parameters 
from the tree amplitudes, while avoiding expanding in 

We adopt a polynomial model for the z-dependence by 



C? M (z) 



(111) 



z{Af M + Af M {6z - 3) + Af M (l 



5z + 5z 2 )}. 



One of these parameters is then eliminated by constrain- 
ing (f M = J dz (f M (z) to its central value obtained from 
experiment. We then eliminate Af* in favor of Cj n to 
obtain 



C BM { 



X) 



2xCf M 



A* M (4x-6x 2 ) 



+ -Af M (x - 6x 2 + 6x 3 ) 
6 



(112) 



Note that the remaining Af M terms in Eq. ( |112| ) must 
integrate to zero. As we will see in the next section, this 
considerably reduces the uncertainty generated by these 
form factor parameters. For M = tt we will simply set 
Af* = 0, thus leaving Af* as t he re maining parameter. 3 
The polynomial form in Eq. ( Ill ) could be justified by 
the SCETn factorization the orem in Eq. (Ha), where it 
is inherited from that in Eq. (10S) at lowest order. How- 
ever we do not view our model in this context, and thus 
do not fix the coefficients Af 1 to values determined by 



Instead we consider Eq. (112) as a model specified 



in SCETi without reference to SCETn-factorization, and 
take Af 1 as parameters to be varied in a suitably large 
range. This ensures that our model for C,f M covers a 
wider range of z-dependence than the restrictive approx- 
imation in Eq. (|4q) would. 



Similar to (f 



(z) we write a model for (f M {z) as 



C^(x) = 2x<£* -Af*{Ax-Qx 2 ) 

+ ^A%(x-6x 2 + 6x 3 ), (113) 

We have taken Cf M (0) = due to the constraint on this 



function derived in section VIII 

For M = p, the simple polynomial model of Eq. (109) 
does not support the value of (x~ 1 ) p = 2.2^0 2 obtained 
from data in section O, unless we include higher or- 
der polynomial terms in Gegenbaucr expansion. Values 
of (x -1 )^/ close to 2.0 require 4> M (x) to peak around 
x = 1/2 with smaller widths. Therefore we choose the 
following model for </> p (x), which has all the desired prop- 
erties 



(j> p (x) = N(a p ) x(l - x) sech 



x - 1/2 



(114) 



Here a p is a parameter whose value is motivated by the 
inverse moment (x~ 1 ) p determined from data in sec- 
tion ||, an d N(a p ) is chosen to normalize 4> p (x) to 1. 



3 The choice Af* = can be justified by isospin and tree-level 
SCETn factorization, but we instead view this choice as part of 
the model. The value of Af K then parameterizes SU(3) viola- 
tion. 
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For ( Bp wc will use a polynomial model like Eq. (112) 



Cj p (z) = 2xCj p + jU* p (x - 6a; 2 + 6a; 3 ) 



(115) 



where for simplicity we take = 0. Alternatively we 



could have based our model for Cj P {x) on Eqs. p4| ) 
and (Q) where it would inherit features of the sech func- 
tion, however we find that using this alternative func- 
tional form does not significantly change our error analy- 
sis. Therefore we adhere to the simple polynomial model 
of Eq. (115). Numerical estimates for the model param- 
eters introduced in this section are presented in the next 
section. 



XII. INPUT PARAMETERS 

Several well determined parameters that are needed for 
our analysis include Q m£ s = 4.7 GeV, mf,(4.7GeV) = 
4.1 GcV, m) s = 1.4 GeV, a s (m b ) = 0.22, fx n (m b ) = 
2.5 GeV, and fiK(m b ) = 2.8 GcV. Defining /m = 
f M /(l GeV) we take f n = 0.131, f K = 0.160, f p = 0.209, 
and from recent lattice data |55| fs = 0.22. We also re- 
quire the Wilson coefficients of the weak effective Hamil- 
tonian, which are known at NLL order |3?|. In the 
NDR scheme taking a s (m z ) = 0.118, m t = 170.9 GeV, 
and m b = 4.7 GeV gives C-i 1 {m b ) = —.316, Cs g {m b ) = 
-0.149 and the NLL results 

Ci-io(m 6 ) = {1.080 , -.179 , .012 , -.033 , .0096 , -.040 , 
4.2 x 10~ 4 , 4.2 x 10~ 4 , -9.7x 10~ 3 , 1.9 x 10~ 3 } . (116) 

In varying /i to estimate uncertainties we will also need 

Ci_i (2m 6 ) = {1.04 , -.104 , .0080 , -.023 , .0074 , -.026 , 
4.2 x 10~ 4 , 2.8 x 10~ 4 , -9.3 x 10~ 3 , 1.3 x 10~ 3 } , 
C 77 (2m 6 ) = -0.281 , C 8G (2m b ) = -0.135 , 

Ci_i (m fc /2) = {1.13 , -.279 , .019 , -.047 , .012 , -.061 , 
5.7x 10" 4 , 7.1 x 10~ 4 , -lO.Ox 10~ 3 , 2.8x 10 -3 } , 
C 77 (m 6 /2) = -0.358 , C% G {m b /2) = -0.166 . 

(117) 

With 2-loop running the MS mass rn b (2m b ) = 3.7 GeV, 
m b (m b /2) = 4.7 GeV, and the chiral-enhanccment pa- 
rameters fi. n {2mb) = 2.8 GeV, fiT T (m b /2) = 2.2 GeV, 
fx K (2m b ) = 3.1 GeV, and fj, K (m b /2) = 2.5 GcV. 



The 7*7 — > 7T° data constrains the inverse pion mo- 
ment, and based on the analysis in Ref. [RTf gives 



-0.03 ±0.14. 



(118) 



For the other linear combination we take — aj = 0.2 ± 



0.3. In our error analysis we do a Gaussian scan over 
these ranges in order to properly take into account the 
correlation in the individual errors of aj and aj, which 
is large. Based on recent lattice data for moments of the 
7r and K distributions we take a 2 = 0.2 ± 0.2 and 
set af = 0. Here the lattice error on was doubled to 
give an estimate for higher moments. For AI = n isospin 
and charge conjugation imply a\ = a\ pp = 0, while for 



M = K wc use 156 



■k,K 
2 0pp 



ipp 

-0.05 ±0.02. 



1. We also take 



77, K 



For simplicity 
0.1 ± 0.3 and 



0.0 ±0.1. For our model of p (a;) in Eq. ( |114| ) wc 



we take 
n K 

use aP = O.ljjv?. 

In section FA we obtained values for the nonleptonic 
form factors ( B7T , (f v , ( Bp , Cf P i from a fit to nonleptonic 
data for the tree amplitudes. Because the uncertainty in 
these parameters are highly correlated we scan over their 
values by doing a Gaussian scan over the range specified 
by the experimental errors in Eqs. ( |2C| ) and ( p3|) and for 
the form factors and 

(^ _1 )mCj M w ith M = w,p. Since 
data is being used for these normalization parameters this 
does not introduce model uncertainty. The choice of the 
remaining parameters Ai introduces model dependence 
to (j(z). We take A\' p = 0, = (0.25 ± 0.30) and 
A p 2 = -(0.05 ± 0.05). We also will use C, B * = 0.0 ± 0.2 
and set A\ x = 0, and A\ 2 = (0.0 ± 0.5). Note that to 
predict the P K7r amplitudes we do not need values of 
C l BK and Cj K , since only pion form factors appear in the 
B — > Kn amplitudes. The kaon form factors are needed 
for B -> KK. 

Finally wc will need values for the model parameters 



appearing in the annihilation amplitudes in section XII 
The three-body decay constants are taken as 
1Q-' 3 GeV 2 , f 3K ~ 4.5 x 10~ 3 GeV 2 , and h p 



4.5 x 
0.13 GeV 2 



from recent QCD sum rule results |57|, |58|, where / 3p = 
iripfj. For the _B-mcson inverse moment appearing in 

the the three-body annihilation amplitude, A 



(i) 

hcann ' 



we 



take Pb = (2.5 ± 1.0)GeV , where the central value 
is consistent with our value for C, B7r using Eq. (^) and 



(109), and the error takes into account the uncertainty 
from the a s (fii) expansion. For the remaining ingredients 
we simply quote results for the necessary moments at 
/j = m b 
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Pic 


— ( — O.Ux l.Dj X 1U , 


Q7T7T 

P-ic 


= U.DoztU.oz , 


PAc 


n i m n no 






/37T7T 

Phcl 


— — 1 .OZ ± U.4Z , 


/D7T7T 

Phcl 


n i o _i_ n in 
= U.io ± U.iz , 


/37T7T 

Phc3 


( O A -L- O 0\ \s 
— ( — IA± L.l) X 


1 n — 3 
10 , 


PhcA 


8** 


= 0.0 ± 5.1 , 


8*1 
Px% 


= 0.0 ± 4.7 , 


8** 

P X 5 


= 0.0 ±0.067, 




8*1 

P x e 


o-nK 
Pic 
















Phcl 


= -1.37 ±0.44, 


ottK 
Phcl 


= 0.13 ±0.12, 


8* K 

Phc3 


= (-2.3 ±2.3) x 


10~ 3 , 


PhcA 




= 0.0 ±6.4, 


o-kK 
P X 2 


= 0.0 ±5.8, 




= 0.0 ±0.085, 




R* K 
P X 6 




= (5.1± 4 ; 2 )xHn 3 , 


8 PP 
PZc 


- -0 11+ 09 

— u.ii_ 03 , 


pz 


= (2.5^-1) x 10- 


-2 

i 




PZi 


= (-3.9t^)xl0- 2 , 


oPP 
Phc2 


= (-1.7i?; 9 6 )xl0- 3 , 


aPP 
Phc3 


= (-1.6 ±1.2) x 


io- 4 , 


aPP 
PhcA 



(119) 

(-4.9 ±1.6) x 10~ 2 , 
0.0 ±0.084, 

(-4.9 ±1.5) x 10" 2 , 
0.0 ±0.10, 



= (i-5±i:f) x io 



1-3 



The values are computed as in Rcfs. E9] and with 
inputs for Ci, fx w , and fix consistent with those given 
above. For the case we used (j>3 P {x, x) = 360xi(l — 

x — x) 2 (x — x)w^, where = —0.20 ± 0.15 is taken from 
QCD sum-rules with an inflated error to account for 
higher Gcgenbauer terms (the relation between our nota- 
tion and theirs is (f>3 P = — $^ p /2). N ote t hat our central 
value of 0.0 for the 8 X i terms in Eq. (HE) indicates that 



over the distribution functions 



we do not have information on the sign of these terms. 
Results for the /3's at fi — mi,/2 and fi — 2rrib are quoted 
in appendix 



p< = 

7T7T 


-!M + 




po = 

7T7T 


-u(( b i - 


^bt)(f*) n 


pc _ 






po _ 


-f K m 




P< = 
PP 


-f P (c c i + 




po = 
p 


-f P ((b c i ~ 


^4)Cf) p , 



(121) 



where /„ = / OT /(lGeV) and f K = / x /(lGeV). In an 
analogous fashion we can define moments for the chiraly 
enhanced penguin amplitudes. For B — ► irir and B — > 
Kn we obtain from Eq. @ and Table | 



XIII. NUMERICAL ANALYSIS 



In this section we make predictions for the penguin 
amplitudes Pm 1 m 2 in the standard model, focusing on 
7r + 7r~, K + ir~, and p + p~ final states. Our sign conven- 
tion for the penguin amplitudes was given in Eqs. ([!]) and 
(^|) . To facilitate comparing the size of various contribu- 
tions we introduce the notation 



MM' 



+ Per + ( 



I pO 

' MM' 

pLann 
r MM' 



"ft 



xC 

MM' 



pxO 
MM' 



r>Gann 
r MM' 



P 



MM' 



bxCx \ 
MM'} 

(120) 



pxC 

7T7T 

pxO 

7TTT 

pxCx 



pxC 



pxO 



P. 



K 7. 



fir Hit f-Brx I XC 

3m B ^ \ ^f«) 



2(9/3) 



3m b 
6m b 



UXC ,uxc 
°3( 9 / 9 ) +0 4( 9 / q ) 



°l(qfq) +0 l(ufu) 



/■Bit 



3itib 
IkPk 
3ms 

fK Pit 

6m,B 



l(qfq) ^ ^2(qfq) 



P)> 
K 



°3(qfq) 



h xc 

°Hqfq) 



h xc 



/■Bi 



A" 



-b xc 

°l(u/u) 



K 



(122) 



The two terms in the first parentheses correspond to the 
leading power terms in line 1 of Eq. (pi]), the second 
parentheses to the chiraly enhanced terms in line 2, and 
P C c corresponds to the long-distance charm penguin in 
line 3. In the last parentheses the first two are LO anni- 
hilation terms from local annihilation and hard-collincar 
annihilation respectively, while the term P^}^ stands 
for chiral-enhanced annihilation. 

The leading power terms can be written as moments 



In Eq. ( |121| ) we have decomposed the leading Wilson co- 
efficients into terms proportional to the two CKM struc- 
tures, 

b ( P = Wb? + X ( t f M = \ { M + X^bl (123) 

where some coefficients (such as and b\) are purely 
from electroweak penguins pCj. Similarly we split the 



Wilson coefficients 



and &f(i?) f° r the chiraly en- 
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hanced amplitudes in Eqs. (|8 



"i(F) 

h x 

°i(F) 



\(f) (X u i \(/) r XC 



\(/) h xu 
A u i(F) 



(124) 



The moments appearing in Eqs. (121) and (|122|) arc 



(biCj M2 )M 1 
( r .\p p 



<iw Ci(u)4> M {u 



(125) 



o 



^BM 2X pp 



J I Mi 
BM 2 



Ah 



1 ,1 

/ 6i(it, (w)Cj M2 (z) , 
Jo 

i 

du a(u)<f>™(u) , 



L/dz6 4 ( W ,z)^HC,f M2 (z), 
o Jo 

1 ,1 

/ dz&^z)^ 1 ^)^ 12 ^)- 
o Jo 



Generically power counting alone gives P^ ~ P^ 7 , where 
the exact size is modified by numerical coefficients. For 
the chiraly enhanced moments the power counting is 



pxC pxO 

MiM 2 Mi M 2 



}XCx 
M1M2 



since C ~ O ~ Cx- 



The penguin annihilation amplitudes can also be writ- 
ten in terms of moments of distributions. Using the no- 
tation in Refs. 29 and pq| the necessary amplitudes are 



P L 



nGann 

TT7T 



in-, 
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f bPb f^T! fit ( ryixiT i fYKIT ^ 07T7T ^ OITTT \ 

^Phd^ Phc2~ 2 Phc3 ~ 2 ci ) 
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2 ^x6 
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K 
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fB0. 



B 
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2^ c4 J 



Kn 



pLann 

PP 



m% 
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m b 



orrK ff^K 



ottK I aizK 

\2 +2 P * 6 



(/3£ + 2/3£ + /?£ 



P 



Gann 



PP 



jBPBf3pfp ( R pp _ opp _ 1 f>PP , 1 a, 

m B m b \ Phcl Phc2 2 Phc3+ * P > 



p/< 

h i:i 



where the /3-moment parameters were defined above in 
Eq. (107) and numerical values were given in Eq. (119). 



To evaluate the remaining penguin amplitudes in 



Eqs. (121) and (|122|) we use the form of the distributions 



from section XL It is useful to write Eq. (125) as integrals 
over short-distance coefficients, i a and j a p, multiplying 
model parameters a a and Ap: 



o/0 

( Ci }m = 4 1c) + E^ c) ^' 

/h c r BM 2 x _ .(4c) >BM 2 \p -(4c) ,BM 2 Mi 

\0 A i,j /Mi — JoO t,J + l_i J /3" Ot 

(0,(9)^(0,0) 

(h CfBM 2 \ _ Alc) r BM 2 , ,-(lc) xBM 2 Mj 

\°1SJ /Mi ~JoO <»J +2^ J /3a A |8 ' 

(o,^)^(0,0) 



(127) 



where a = 0,1,2,4 with a A J = (1, af, a M , af), and 
p = 0,1,2 with Af M = {(f M ,A? M ,A BM ). This step 
is useful because the short-distance coefficients, i a , j a p 
are integrals which can be evaluated numerically inde- 
pendent of the choice of the model parameters. This 
makes it easier to propagate errors from parameter un- 
certainties into the final amplitude predictions. It also 
makes it possible to study the short-distance uncertain- 
ties (such as the /^-de pendence) directly in terms of i a 
and jf3 a . In Eq. (127) we have separated out the dom- 
inant term from the sum. Since our values of Cj T and 
C B7r are extracted from independent experimental data, 
these dominant terms in the penguin amplitudes become 
model independen t. Fo r the chiral enhanced amplitudes 
the analog of Eq. (127) is 



/„XC \PP 
\ C l(qfq)/M 



/„YC \pp 
\ C 2(qfq)/M 

luxe r BM 2 \PP 
\°3(g/g)^J /A/i 

llXfi sBM 2 \pp 
\°4(qfq)^J I Mi 

luxe /-BM 2 \ 
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$X{ufu)^X 2 ) M i 
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"x/3 fl « 



where /3,7 = 0,1,2, pjf = {l,a^ p ,a^ p } and = 
{Cx A/ , ^i M , ^x2 M }- In tcrms of the J' fe ' ^ coefficients, 
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Eqs. dm]) and djU) are given by 
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and 
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(130) 



where a sum over a = 0, 1, 2, 4 and (3, 7 = 0, 1, 2 is un- 
derstood. 

Evaluating the short-distance integrals at zeroth order 
in a s with the Cj's in Eq. (|116D the i's and j's are 



l ( 4c ) x 10 3 
i( lc )xl0 3 

J (4c) xl0 3 

j (lc) xl0 3 



(-28.4, 0, 0, 0) , 
(-1.96, 0, 0, 0), 

-11.1 17.3 17.3 17.3 



-16.4 -14.5 -14.5 -14.5 





(131) 



Note that the short-distance "i" coefficients for £ are 
comparable in size to the short-distance "j" coefficients 
for (j M (z). For the chiraly enhanced integrals we find 



k^ 
k^ 

die) 



xlO 3 = 
x 10 3 = 



xlO 3 = 



/2c) x 1Q 3 



(3c) 



xlO 3 = 



/4c) x 1Q 3 = 



(218, 0, 218), 

(-3.38, 0, - 3.38), 

-12.3 69.1 -69.1 -69.1 
34.5 -34.5 34.5 34.5 


61.9 -57.9 57.9 57.9 
-29.0 29.0 -29.0 -29.0 


218 218 



-3.38 -3.38 





(132) 



Relative to the size of i^ 4 ^ and j( lc > 4c ) the enhanced size 
of the k^ and the short-distance coefficients is 
quite striking. Comparing the matching coefficients in 
Eqs. (HQ and (|l]f||) we sec that the combinations of 
coefficients from Hyy are similar in size (C34 versus C$fi). 
However, the k^ lc ^ and £^ lc ' moments are enhanced by a 
factor of ~ 6 due to the inverse moment fraction factor 
1 /uu. This numerical factor provides additional enhance- 
ment beyond the numerical enhancement in /jm/wi, and 
is the essential reason why the chiraly enhanced penguin 
amplitudes are numerically important. 

Next we evaluate the short-distance integrals i^ 4c \ 
j^ ic \ and /c( lc ) up to order a s , by including the one- 
loop results for b± and C4 given earlier in section 0. 
Where known we also evaluate the chiral enhanced short- 
distance integrals up to 0(a s ) (from Eq.([32"|)). 



At the scale [i = mj, we find 



i(4c) 



xlO 3 = 



j^xlO 3 
fc (lc) xl0 3 



(-36.2-i9.58,4.49+i8.12,10.5+i2.06,4.42-i2.36), 

-14.0-i4.52 27.6+^4.47 28.1-i0.15 25.8-i0.55 
-1.05-i3.63 -2.30 + i2.10 1.32+il.95 -0.20-il.64 
0.15+i0.04 -0.10-i0.08 0.08-i0.14 -0.07+i0.14 

(281+i54.7,8.6-i69.1, 240+i41.1) . 



(133) 



We will also analyze how stable our results are to variations in /1. For the LO results in Eqs. (131) and ( |132| ) a change 
in /1 simply reflects changes in the C,(/x) and so will not be shown. At NLO in the pcrturbative expansion we find for 
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H = m b /2 = 2.35 GeV 



x 10 3 

/4c) x 1Q 3 



(-38.9-il2.3, 7.46+il0.5, 15.2+i2.65, 7.37-^3.05), 

-12.6-i6.67 34.6 + ^6.12 36.3 + i0.45 33.l-il.23 
-0.92-i4.01 -2.51+i2.30 1.40-M2.14 -0.13-il.82 
0.20 + i0.09 0.17+i0.09 0.08-t0.21 -0.08 + i0.20 

(359 + i70.6, ll.l-z89.2, 306±i53.0) , 



while for (i = 2m& 



fe (ic) x 1Q 3 = 
9.4 GeV wc find 

j(4c) x 1Q 3 = (_32.5-i7.57, 2.81+i6.42, 7.58+il.63, 2.75 - £1.87), 



/4c) x 1Q 3 
/fc (lc) x 10 3 



-13.2-i3.10 22.7±i3.31 22.6-i0.45 21.0-i0.16 
-1.05-i3.19 -2.05+il.85 1.20-M1.72 -0.22-il.44 
O.ll + iO.Ol 0.06 + i0.08 0.07-i0.09 -0.07+i0.10 

(228 + i43.3, 6.8-i54.7, 195 + i32.5) . 



(134) 



(135) 



From Eq. (133) we observe that these a s corrections in- 



duce imaginary contributions which are often appreciable 
since the a s Ci^ terms can compete with C 3 _ 6 . For ex- 
ample, the imaginary part of Jqq C ' determined from our 
result for the one- loop matching given in Eq. ([52]) , is 
~ 30% of the real part. 

Because we have neglected terms a s C3_6 we must also 
neglect the fi dependence of (, £j, and the <^'s for con- 
sistency. These terms induce a a s ln(^t) that multiplies 
the tree-level penguin coefficients involving G^-q and 
are hence compensated by a s ln(/i)C3_g corrections to 
the short-distance coefficients. The dominant coefficients 
have a = (3 — 0. At zcroth order in a s the central val- 



ues for the coefficients i 4c ^ and j^ ' vary by ±30-50% 
when we take fi — mb/2 and fi — 2m;,. We find this 
change is reduced to < 10% at NLO. At LO the chiraly 
enhanced fcg lc ^ varies by ±35-55%, and this is reduced by 
about a factor of two, to ±20-25% at NLO. The imagi- 
nary parts first appear at 0(a s (/j,)), and exhibit a ±20- 
30% range for and ky lc \ and ±30-50% dependence 
for j'oq^. The LO coefficients ^ 00 c ' 3c,4c ^ also have a size- 
able /x-dependence (20-50%) and it will be important to 
compute their a s corrections in the future. Below we will 
take this residual scale uncertainty as a way of estimating 
the size of missing higher order perturbative corrections 
on our final result. 



(4c) 



On the other hand the electroweak coefficients 



.(lc) 



and £^ have only ~ 3% /i-dependence at LO, consistent 
with our expectations that the NLO corrections to this 
term are small. This reflects the fact that the correspond- 
ing \i dependence occurs in a NLO penguin diagram with 
photon exchange, whereas the leading order Wilson co- 
efficients are generated by both photon and the larger Z 
exchange. The corrections to the electroweak coefficient 
»q at LO is larger as a percent (40-50%), however the 

.(lc) 

i\ is tiny to begin with, since at LO its proportional to 



the numerically small combination Cm + Cg/N c . Thus we 
do not expect our neglect of these one-loop electroweak 
corrections to have a large effect . 

In table || wc pre sent numbers for the penguin am- 
plitudes in Eq. (12C), showing separately the tree- level 
and a s corrections. The errors shown in the table in- 
clude only input parameter uncertainty, and are com- 
puted with Gaussian scans for the errors in the model 
parameters given in section XII. Despite having a num- 
ber of hadronic parameters, we observe a relatively small 
model parameter dependence in the first four columns for 
7T7T and Kir. This occurs because there is only a small 
dependence of the penguin amplitudes on the shape of 
Qj M {z) and 4> AI (u). The normalization terms, ig 4 ^ and 

(lc Ac] 

j 00 : , give the dominant contribution to the amplitudes 
and the corresponding model parameters were fit to inde- 
pendent data to reduce their uncertainty. Furthermore, 
at LO the only shape parameter dependence comes from 
and for the pion this parameter is quite well 
known (which in our error analysis is accounted for by 
taking into account an important correlation in a\ and 
of). For example, we find that further doubling the error 
bars on the shape parameters Afl x 2 only effects the last 
quoted digit of the error bars on the penguin amplitudes 
m table H 

Since the second and third rows of jj^ are 

small, the dependence on the shape of (f M (z) is very 

small. Though the coefficients jg 4 *^ , are similar in 

size to ioo , their contributions are suppressed by the 
small af^. One entry is very sensitive to the , namely 

the a s Ci t 2,8g contributions to Pj^. Here the dominant 
term is 



io 3 p£;; 



] K C,T [( 2 - 9 + * 4 - 5 ) - ( 10 - 3 + i4.47)af 



- (10.8-i0.15)af] , (136) 
and exhibits a large cancellation in the real part for the 
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Penguin Amplitudes, Pm 1 m 2 x 10 4 






Penguin Annihilation, Pm 1 m 2 x 10 4 




P c 

J 7T7T 


po 

7T7T 


pxC 

J TTTT 


pxCj 

J 7T7T 


pxC x 

* TTTT 




pLann 

7T7T 


nGann 

7T7T 


pLannx 

J 7T7T 


C3-10 


3.58 ±1.02 


3.34 ±0.88 


4.41±1.78 


4.51 ± 1.71 


0.00±1.03 










Q!sCl,2,8s 


(0.86±0.25) 
+i(1.08±0.32) 


(0.32±0.26) 
±i(0.53±0.13) 


(1.21±0.37) 
+i(1.10±0.40) 






a s Ci_io 


-1.46±0.88 


0.15±0.08 


0.00±5.00 


P^XIO 4 : 


P c 


pc, 

Kir 


pxC 
r Kir 


pxC./ 

r Kir 


pxCx 

r Kit 




pLann 
^ Kir 


pGann 
-'Kir 


pLannx 


C3-10 


4.37 ±1.25 


4.00 ±1.02 


6.02±2.42 


6.15 ±2.33 


0.00±1.18 










&sCl,2,8g 


(0.86±0.40) 
+«(1.38±0.40) 


(0.11±0.35) 
±i(0.70±0.18) 


(1.66±0.50) 
±i(1.50±0.57) 






a s Ci_io 


0.26±0.14 


0.20±0.11 


0.00±8.03 


p; p xio 4 : 


P c 
pp 












nLann 


pGann 
PP 




C3-10 


14.8 ±3.5 


1.64 ±2.99 


















(5.59t?;°f) 


(0.36±0.66) 
±i(0.28±0.49) 








Q s Ci-io 


0.65l° ; 6 2 ? 


fj rjrj+0.22 





TABLE II: Numerical predictions for the penguin amplitudes, Pmm, from the factorization theorem. The results were split 
into terms generated at 0{oP s ) and 0(a s ) in the short-distance matching coefficients. In each row the theoretical predictions are 
broken down by contributions from the C, BM and Cj M terms in the leading power amplitude, terms C, BM , C, BM , C, BM from the 
chiraly enhanced part of the amplitude, and terms from standard penguin annihilation, three-parton annihilation, and chiraly 
enhanced annihilation. The errors shown are uncertainties propagated from input parameters as described in the text. 



value = 0.2 (explaining its large percent uncertainty 
for = 0.2 ±0.2). Some cancellation is also evident for 
P%£- However, overall these are both small contributions 
to their respective penguin amplitudes. This type of pa- 
rameter dependence does not appear in other terms, and 
we find that it does not significantly effect the final result. 
Thus even though our model parameters vary over a large 
range we have fairly robust central values for individual 
contributions to the leading order penguin amplitudes in 
table y. The uncertainty in the 3rd and 4th columns for 
the chiraly enhanced amplitudes is also reduced by our 
knowledge of the normalization of £ and £j, and is a bit 
bigger than the first two columns due to the added un- 
certainty from (f>pp(u). The 5th column involves the new 
form factor ( BM (z), where we do not have information 
about the sign, and hence zero central values. 

In the 7T7T and K n entries in table || we also observe 
that the contributions from Q BM and (f M arc similar 
in size. This is a reflection of the fact that there coeffi- 
cients are similar numerically, and is in agreement with 
the power counting £ BM ~ (f M . In determining the 
errors associated with these parameters it was quite im- 
portant to take into account the correlations, as already 
described in section [n]. Also, as mentioned above, the chi- 
raly enhanced penguin amplitudes compete numerically 
with the leading power amplitudes due to the presence of 
the enhancement by the 1 juu momentum fraction factor 
which generates a numerical factor of six. For example, 
we have 



fi(j _i_ £>xC 

TTTT * TTTT 



~/ w Cr(28 + 215 
C3-10 ^ 6m 



(137) 



where the large numerical value 215 is generated by this 
enhancement. 



Examining the annihilation amplitudes we see that 
P^ nn + P^ nn is suppressed by a factor of 1/5 relative to 
P^n + P£%, and so is of the expected size for this power 
correction, namely 



pLann 

TTTT 



pG 



ann 

7T7T 



A 



P c 

1 7T7T 



pc, 

1 7T7I 



_ a s (m b ) 



QCD 



(138) 



The same conclusions hold for Kir. On the other 
hand the chiral enhanced annihilation terms P^ nnx and 
P K ^ nx have much larger parameter uncertainty, and we 
are not able to draw definite conclusions about the size 
of these terms. In fact they provide the dominant pa- 
rameter uncertainty for the -k-k and Kn channels. 

For B — > pp decays, our analysis was slightly differ- 
ent from the pseudoscalars as we used a non-polynomial 
model for <fi p (x). Here the errors are dominated by the 
uncertainty in ( Bp , Cj P : an d a-p- The uncertainty from 
the shape parameters Ai are negligible in comparison. 
Since current data prefers a central value for £j 9 sig- 
nificantly smaller than that for C, Bp this same hierarchy 
is observed in the penguin amplitudes. The size of P£ p 
is enhanced by f p and the p-form factor in comparison 
to P^tt- Due to the absence of chiraly enhanced contri- 
butions the dominant parameter uncertainty comes from 
experimental uncertainties that propagate into the errors 
for C Bp andCj p . 

the individual contributions 



In tabic III 



we sum up 

from the leading power, chiraly enhanced, and annihila- 
tion penguin amplitudes, to obtain P LO , P x , and P ann 
respectively. To perform these sums we do separate 
Gaussian scans for the total penguin amplitude since this 
provides the simplest way of propagating correlated pa- 
rameter uncertainties. This also explains why the central 
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P LO x 10 4 


P x x 10 4 


pann ^ 1Q 4 


ptotal x 1Q 4 


Ar P £*io 4 

( 7 = 59°) 


^xlO 4 
(7 = 74°) 


P^ pt x 10 4 
(7 = 59°-74°) 


B — > 7T7T 


(8.10±0.63) 


(10.2±2.9) 


-1.31 ±5.08 


(16.9 ±5.9) 


(18±9) 


(44±6) 




±£(1.61±0.21) 


+i(1.10±0.39) 




+£(2.71 ±0.45) 


-£(29±6) 


-£(29±6) 




B — > Kit 


(9.34 ± l.OO) 


(13.8 ±3.9) 


0.46 ±8.03 


(23.6 ± 9.0) 






±(48±4± 10) 


±£(2.08 ±0.25) 


+£(1.49 ± 0.57) 




+£(3.57 ±0.62) 






-£(22±7±4) 


B^pp 






0.871S;S 




-(29 ±26) 


(38 ± 23) 




+i5.68tlf 7 






+£5.68±?; 4 | 


-£(8 ±18) 


-£(8 ±18) 





TABLE III: Numerical predictions for the short-distance penguin amplitudes at leading power, P LO , from chiraly enhanced 
terms P x , and from the annihilation amplitudes in Refs. pq, E9|. The sum of these contributions p* *^ is the total short- 
distance result from the factorization theorems discussed in the text (long-distance terms are discussed in the text). The last 
three columns show current experimental data. Comparing them with p total shows an order of magnitude short-fall for the 
imaginary part. 



values arc not precisely the mean from table O, due to 
small non-linearity effects in the parameter dependences. 
The correlation in input parameter uncertainties must 
be taken into account to get the errors shown here. The 
three amplitudes in the first three columns of tabic III 
are then added together to get the total theoretical con- 
tribution, p total . These total values can be compared to 
the experimental values in the last three columns. The 
uncertainty shown only includes the variation of param- 
eters from the Gaussian scans. For the first column the 
displayed errors are dominated by the uncertainties in 



for all theoretical uncertainties we find 



u 2 -r u 4 , af 2 i ( B7r 7 Cj^i an d for B — > pp those in a p , 
C Bp , and (j p . The effect of other parameter uncertain- 
ties is quite small. Even the dominant uncertainties are 
small due to our proper account of parameter correla- 
tions and use of experimental data. Also due to our fit 
procedure the errors from £ and £j will decrease with im- 
proved measurements of the tree amplitudes (which come 
from improved branching ratios and CP-asymmetries). 
In p total the uncertainty from the parameters in the chi- 
ral enhanced annihilation by far dominate the errors for 
B — ► 7T7T and B — > Kir. 



In addition we can estimate the uncertainty from deter- 
mining the hard coefficients by varying p G [rrib/2, 2mb] 
For the real parts this gives an additional 



total 

7T7T 



+ 15% 
-12% 



+7% 
-9% 



° o uncertainty for , and 



uncer- 

+9% 



-10% 



-19% 

for P^l a \ and ^f 2 % uncertainty for P\ 



total 
PP ' 



tainty for P\ 

uncertainty for Pp° p tal - For the imaginary parts we find an 
additional +^ uncertainty for -P*° tal , uncertainty 

Finally we as- 
sign a generic 20% uncertainty to the final p total results 
to account for the fact that we have given only a par- 
tial treatment of l/m& corrections, but do not foresee a 
reason why the untreated corrections should be enhanced 
over the power counting estimate. Thus with an estimate 



P 



total 



p, 



.total 
K-k 



p, 



total 
PP 



(16.9±5.9l}?±2.0±3.4) 

+ i(2.71±0.38±;^±0.33±0.54) , 
(23.6±9.0^j±2.8±4.7) 

+ i(3.57±0.53±;^±0.43±0.71) , 
(23.3±^;^±2;3±2.8±4.7) 

+ i(5.68^;^ij;^±0.68±1.14) . 



(139) 



The first errors are from input parameters and are dom- 
inated by chiral-enhanced annihilation for B — > 7T7t, Kit. 
The second errors are our estimates of higher order per- 
turbative corrections (the /x- variation) . The third terms 
are errors from \V u b\ which propagate through the form 
factors and hence can be added as a ±12% uncertainty. 4 
Finally the fourth errors are a generic 20% that we add 
for unknown power corrections. 

For 7T7T the real part of the amplitude in Eq. (139) 
agrees with the data in table III for 7 = 59°. However, 
the same is not true for Kir, nor even for -kit if 7 = 
74° (which is the value preferred by SU(3) and SCET 
power counting which predicts P W7r ~ Pk-k |Tq])' Here 
the disagreement with data in the real part is at the level 
of factor of two. 

On the other hand the imaginary part of the short- 
distance prediction for F"" and P Kn are much smaller 
than the corresponding experimental values and have the 
opposite sign. Due to a numerical enhancement P^j lM 

and PmiM 2 arc °f same Slze as t ue leading power contri- 
butions to the amplitude, but as wc have demonstrated 
by deriving an SCETi factorization theorem, these terms 
are real at zeroth order in a s . After taking into account 
all theoretical uncertainties in our analysis, we conclude 
that it is not possible to match the P imaginary parts 



4 We have increased the 7% error on \V u b\ quoted by HFAG 
which we consider to be overly optimistic. 
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obtained from experimental data. Therefore the large 
phase of the penguin relative to tree amplitudes can only 
be explained by long distance charm contribution, P cg , 
within the standard model, or by contributions from new 
physics. 

If the remainder is generated by long distance charm 
contributions, then we can determine what values of 
Pec 1 2 reproduce the experimental data. This gives 

p™ = { (l±H)-i(32±6) (7 = 59°) 
1 (27 ±8) -i(32±6) (7 = 74°) ' 

P£ 77 = (24 ± 14) - i{26 ± 8) , (140) 

where we have added the experimental and theoretical 
errors in quadrature. Thus a long-distance charm pen- 
guin with substantial imaginary amplitude is one pos- 
sibility for reproducing the data. This explanation was 
favored in Refs. || [2^, ^3), and the analysis here makes 
the required size of these long-distance terms fairly pre- 
cise. In the next section we contrast this long distance 
standard model explanation with the more exciting pos- 
sibility of a new physics contribution. An additional 
test of the penguin amplitudes can be made from study- 
ing the channels B~ -> K°K~, P° -> K°K°, and 
B° — > K + K~ which get contributions from penguin and 
annihilation/exchange type diagrams. Branching ratios 
for these channels are available |nj. Since for B — > KK 
we do not have enough experimental information to fix 
Q BR and Cj R we resort to SU(3). We can apply SU(3) 
directly at the level of SCET amplitudes as discussed in 
Ref. Jl§, and it implies that P KK ~ P K7r up to the 
small penguin annihilation terms. Using the experimen- 
tal value of P^iKn) from table H this value of P KK 
reproduces the data for Br(B~ -^K®K~) = 1.36 ±0.28 
and Br(B° K°K°) = 0.96 ± 0.20 from HFAG §. 
The channel B° — > K + K~ does not get contributions 
from the penguin amplitudes P^ , P^ J , nor the chiraly en- 
hanced penguin amplitudes as can be seen from Table |. 
It does get contributions from annihilation, but not from 
the potentially sizeable chiraly enhanced penguin anni- 
hilation, p Lm ™x i as is c i ear f rom Table IV of Ref. [§9). 
Hence the small observed value Br(B° — > K + K~) = 
0.15 ± 0.10 is consistent with the size of the annihila- 
tion results for p Lann an d F "™ in table || The size of 
these amplitudes is also consistent with the power count- 
ing estimate of A/ nib suppression relative to leading order 
terms. 

The experimental errors in P pp are too large at this 
time to draw strong conclusions, but it is interesting 
to note that the positive sign for the real part of the 
short-distance standard model penguin prefers values of 
7 larger than 59°. 

Our numerical results for the penguins can also be com- 
pared with earlier analyses in the BBNS (l?], [27j] and 
KLS |25|, ^6) approaches where light-cone sum-rules are 
used for the hadronic parameters. The BBNS analysis 
also gives numbers where P x ~ P LO , and gives small 
short-distance imaginary parts in P L °. However, indi- 



vidual central values differ from ours due to their different 
method for dealing with input parameters and their use 
of an expansion in a s (/j,i) at the intermediate scale for the 
LO penguin and chiral enhanced penguin contributions. 
Also a larger (complex) range of annihilation amplitudes 
was adopted in Ref. J27J , with a non-perturbative strong 
phase that can be chosen to fit the data. In the KLS 
approach it is more difficult to compare individual con- 
tributions, but generically the penguin amplitudes are 
somewhat larger, and have a large strong phase from an- 
nihilation graphs. The most prominent feature in both 
comparisons is that our parameter errors in P LO and P* 
are significantly smaller than earlier results, due to our 
use of tree amplitude data to determine the hadronic pa- 
rameters. From our numerical analysis of annihilation 
amplitudes, together with power counting arguments it 
appears that nonperturbative charm loops are the most 
likely culprit for a missing long-distance contribution to 
the amplitude. 



XIV. PENGUINS CONTRIBUTIONS FROM 
NEW PHYSICS 

There has been a lot of discussion about the possibil- 
ity of new physics in nonlcptonic P-decays (for exam- 
pie §jj, H II 0, 0, |§ g, [r§ [§, 0, @). The 
precision achieved for the computation of the standard 
model penguin amplitudes in tables [n] and [n|, and their 
lack of concordance with the experimental results, make 
it interesting to reexamine the role new physics contribu- 
tions may play. In this section we aim to look at general 
features the new physics contributions should have, and 
do not attempted to explore this topic in specific models. 

Lets consider adding new physics contributions to 
the nonlcptonic amplitudes. Since new flavor-changing 
physics is likely to be heavy we can suppose that upon 
integrating out the short-distance new particles we gen- 
erate a set of operators whose amplitude is parameterized 
by a CP-even matrix element N and a CP-violating phase 
</>, A NP = Ne 1 ^. Here N contains the strong rescattering 
phase for the amplitude, and e 1 ^ has CP violation that 
need not follow the CKM paradigm. 

In order to fit the data, i.e. contribute to Im(P), we will 
demonstrate that Ne 1 ^ must have a non-zero CP-even 
strong phase. Given this we may ask whether a small 
strong phase in N can be enhanced by a large new source 
of CP violation, or by some other new physics effect. Wc 
will see that there is a strict bound that prevents us from 
enhancing Im(P) without having large Im(N). 

To study these points we follow Ref. and use the 
fact that we can decompose any new physics amplitude 
into terms that simply shift the CP-even standard model 
amplitudes in Eq. (]]]). For example, we can decompose 
any Ne 1 ^ to make it look like terms appearing in B — > 
7T7T, Kir, pp: 

Ne^ = Ni + N 2 e'^ , (141) 
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where the first terms acts like the Ai d ' s ' term and the 
second like X^f' s \ Here only <fi and the standard model 
phase 7 change sign under CP, while Ni = + iN[ 
and N 2 = N^ + iN^ are CP-even. Adding Ne^ terms to 
the SM amplitudes, we see t hat N 12 simply shift the SM 
amplitude parameters. Eq. (141) was used in Ref. [jF7| to 
point out that it is not possible to observe new physics 
in penguin amplitudes in decays like B — > 7T7t, Kit with- 
out having information about the SM penguins that goes 
beyond isospin. Given the computations of the SMpen- 



guins in the previous section, we can use Eq. (141) to 
explore how new physics effects can appear. To gener- 
ate large Im(P) in our phase convention we need large 
Im(iVi) and/or large In^A^). 

Being CP-even the parameters A1.2 act like strong 
interaction amplitudes, despite the fact that they con- 
tain sho rt-distance CP-violating parameters. Solving 
Eq. 



141) gives 



N 



N 



Im(JV) 
sin 7 

Re(iV) 
sin 7 



sin(7 + <f>) 



Im(JV) 
sin 7 



sin(0) . 



sin( 7 + 0) , Ni 



Re(A) . 

sin(0J 



SU17 
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Hence the shift to the imaginary part of the standard 
model amplitudes is zero if Im(iV) = 0. Furthermore we 
have the bounds 



\N[\ < 



Im(iV) 



sin 7 



1^2 I < 



lm(N) 



sin 7 



(143) 



The SM value of 7 is not small (sin 7 ~ 0.9), so these 
bounds imply that enhancement in N[ or N 2 requires 
large Im(7V), and hence a large strong phase for this new 
physics amplitude. Thus given sin 7, no enhancement of 
the effective strong phase can occur due simply to new 
sources of CP-violation. This conclusion does not appear 
to be changed if one or more new physics amplitudes are 
added in the various standard model decay channels. 

The CP-even phase in N will be generated by strong 
rescattering, and it is useful to consider ./V as an ampli- 
tude generated by new dimension-6 four-quark operators 
not present in the SM. Our analysis of SM four-quark 
operators gave power suppressed non-perturbative strong 
phases and small strong phases from hard penguin loops, 
so we might speculate that the same would be true for 
four-quark operators with non-SM symmetry properties. 
In th is ca se the imaginary part of N will be small, and 
Eq. (143) implies that adding new physics will not sig- 
nificantly improve the situation with Im(P). One might 
think that the inclusion of new physics into the process 
of extracting a value for P could mollify the need for a 
large imaginary part. However, a simple analysis, say in 
the 7T7T modes, shows that the existence of an N with a 
small imaginary part can not lead to penguin completion. 
It will simply shift the meaning of the real parts of the 
tree and penguin amplitudes in the fit, with only a small 
change to the meaning of the imaginary parts. 



Thus for new physics to play a significant role in the 
observed Im(P) we need to find a large imaginary part 
for an N from analyzing an operator not generated by 
the standard model. Though there is no reason to ex- 
pect an enhanced short-distance contribution, this is a 
logical possibility which deserves further study. Signifi- 
cant new physics contributions could in fact be obtained 
by modifying the coefficient of four-quark operators with 
charm quarks, since then a large long distance charming 
penguin amplitude could provide the necessary contri- 
bution in Im(iV). It might be interesting to attempt to 
construct explicit new physics models of this type which 
are not ruled out by other constraints on flavor changing 
neutral currents. Thus it seems to be quite a challenge to 
complete the penguin without the aid of a long-distance 
contribution. 



XV. DISCUSSION AND CONCLUSION 

Let us now address the question raised in the title. 
The results in table III show a lack of concordance be- 
tween the theoretical prediction for the short-distance 
standard model penguin contributions and the extracted 
value for the penguin amplitude. 5 Chiraly enhanced op- 
erators substantially increase the penguin contributions, 
but they are not able to generate the necessary imagi- 
nary pieces. Thus it would seem that the shortfall must 
be due to cither the long distance charm or new physics. 

Before addressing these possibilities we must be as- 
sured that the assumptions leading to this conclusion are 
justified. Our theory predictions for SM penguins assume 
that the expansion in powers of A/mt, is trustworthy, 
since the convergence of this series is a necessary criteria 
for factorization to apply. The experimental extraction 
of P n7r and P pp relies on isospin, and hence is quite ro- 
bust. The penguin extraction for the Kir system relies 
slightly more on the factorization (the A/mb expansion) 
since we use factorization for the tree amplitude T K ™ . 

What evidence do we have that the large mass ex- 
pansion is indeed converging? The factorization theorem 
for color allowed B — > D^*>M~ decays (proven in Jf9|) 
agrees with data with the expected accuracy. For color 
suppressed charmed decays, B -> D^M° the SCET 
prediction for the strong phases |3J 



is in good ac- 
cord with the data for many channels, which provides a 
non-trivial test of the large energy expansion. The same 
expansion is also used in analyzing the photon cut de- 
pendence of B — > Jf s 7 |Hj and for the analysis of \V u b\ 
from B — ► X u lu |p2| , where power corrections appear 
with the expected size. One might object that these last 
two examples are inclusive, summing over states up to 
/i 2 ~ TO{,A. However, our analysis is quite similar, since 



5 As shown in Eq. for this conclusion uncertainties in the 

weak phase 7 are irrelevant for lm(P n T r ), but not for Re(P 7r7r ). 
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the factorization theorems we use do not attempt to fac- 
torize physics below /i 2 ~ to^A, and instead retain it as 
form factors. Due to experimental cuts an analysis of 
B — > X s £ + £~ data will also rely on this type of expan- 
sion f§. 

More direct evidence for our methodology for analyz- 
ing charmless nonleptonic decays comes from successes 
in the exclusive modes themselves. In [^5[ [l||], a com- 
plete list of the predictions for branching ratios and CP 
asymmetries was given, by using the data to fit the 
unknown hadronic parameters (including long-distance 
charm penguin amplitudes). The theory fits the data 
quite well, with all of the theory points falling within 
1-2 a of the data. (The only significant exception is 
the ratio of A C p{K ~n )/A C p{K ~7r+) where the sign 
disagrees with the data.) It is interesting to note that 
SCET predicts certain asymmetries to be negative while 
the current experimental central values are positive. The 
factorization theorem for charmless nonleptonic decays 
also gives a prediction for |V^i,|/+(0) given in Eq. (20), 
which is in good agreement with the recent extractions 
based upon dispersion relations |84|, |85| utilizing lattice 
data |6^, |63|, ^3j|. (Using Hill's 8 parameter Q an 
analogous test will be possible for £j ^ with future exper- 
imental improvements on the B — > tt£D spectrum.) Note 
that for all of these successes the penguins were fit to the 
data and any deviation from the short distance predic- 
tion was absorbed into the long distance charm piece A cc . 
Thus these successes do not directly imply convergence of 
the large energy expansion for the penguin amplitudes. 
However, from the point of view of QCD there is not 
much distinction between short distance tree and penguin 
contributions. Although the pattern of contributions to 
each of their amplitudes differs, the A/rrib expansion for 
each type of contribution involves very similar hadronic 
physics. An exception occurs for charm quarks, where 
the non-relativistic region and poorer convergence of the 
A/m c expansion may play a role. 

For the penguin amplitudes we can see from table III, 
that the chiraly enhanced power correction is of the same 
order as the leading order penguin contribution. It is in- 
teresting to understand the origin of the enhancement 
for these power corrections. First off, the chiral conden- 
sate gives an enhancement of a factor of three ^7|. As 
we discussed the chiraly enhanced contribution also has 
a Wilson coefficient which gives an added numerical en- 
hancement compared to leading order penguins by ~ 6, 
coming from a factor of 1/ (uu) . One should worry that 
there could be higher dimensional operators which are 
chiraly enhanced as well. However, for these operators to 
be as large as the leading chiraly enhanced contribution 
they would have to have additional enhancement from 
their coefficient function. At present no such sublcading 
operators are known to exist, but further investigation is 
warranted. 

From the power counting and form of the factorization 
formulas for the nonleptonic amplitudes, terms with a 
long distance phase (outside of long-distance charm am- 



plitudes P cc ) arise from contributions which are down by 
a(/Xj)/7r relative to the corrections considered in this pa- 
per, sec |36|. In principle these corrections, which come 
in at order A/ 'm(,a(/Xj) / V , as well as A/xm/w 2 ct(/Xi)/7r 
(the chiraly enhanced pieces), could account for the pen- 
guin deficit. However, this would be in gross violation of 
the power counting. Even if the expansion in a(fii) were 
very poorly behaved, which seems not to be the case in 
the calculations performed to date [^0), these contribu- 
tions could still not make up the deficit, as they would 
be expected to be the same size as the chiraly enhanced 
penguin annihilation (at best), shown in the last column 
in table II. Of course if the chiraly enhanced annihilation 
were truly as big as the lower order terms in the power 
expansion, which our error in the table allows, then we 
would question the whole power expansion in the pen- 
guin sector. However, to push the penguin annihilation 
to the limits of our errors one needs large deviations from 
naturalness. 

Two possible resolutions are, new physics and long dis- 
tance charm. Let us consider the former possibility first. 
As we have shown in the previous section, introducing a 
large CP violating phase from new physics, does indeed 
have the effect of mimicking a CP conserving imaginary 
penguin, however its size is bounded by the strong phase 
induced by QCD. Thus, given that we have shown that 
such imaginary pieces (modulo A cc ) are small, it would 
seem to be a challenge to complete the penguin using new 
physics. One open possibility is that the new physics gen- 
erates new operators not present in the standard model 
electro-weak Hamiltonian, that generate large imaginary 
parts when matching onto SCETi. However, given our 
experience matching the standard model operators, there 
is no compelling reason to believe that such a scenario is 
likely. 

In addition, generically, the new physics is constrained 
to only arise in certain operators. In particular, we note 
that the new physics would not fall under the rubric 
of Minimally Flavor Violation |6(], pT^ , since there are 
strong constraints on e' in the kaon system. This in itself 
is not a problem as one might expect the new physics to 
couple differently to the third generation, given the top 
quark mass. Furthermore the new physics should leave 
the AB = 2 operators responsible for B — B mixing es- 
sentially unscathed. It would seem to be an interesting 
challenge to build a model which accomplishes these goal 
without fine-tunings. 

Long distance charm contributions are perhaps the 
most compelling expla natio n for the penguin deficit. As 
was shown in section XIV current data appears to re- 



quire a sizeable long-distance strong phase, such as the 
long-distance charm amplitude described in section IX. 



Moreover, the long distance charm has the potential to 
explain another discrepancy with the data ]87[ |, namely 
the deficit of transversely polarized vectors in the de- 
cay (j)K* channel ||. In SCET one does not generate 
any leading operators which produce transversely polar- 
ized vectors. This suppression follows from simple chi- 
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rality arguments pifl . To derive an amplitude factor- 
ization formula for the long distance charming penguins 
which generates transverse polarization was beyond the 
scope of this paper. However, it is simple to see that 
we would expect transverse polarization by noting that 
the helicity arguments mentioned above no longer apply 
because the valence quarks which make up the mesons 
are no longer produced on the light cone. Moreover, a 
post-diction of our SCET analysis method would be that 
one would expect a large transverse polarization fraction 
in the <pK* channel, but not in the pp, simply because 
the latter is tree dominated while the former is penguin, 
and hence A cc dominated. Thus it would seem that the 
long distance charm contribution can explain both the 
penguin dearth as well as the transverse polarization in 
<f)K*. Whereas a new physics scenario would seem to 
need some organizing principle which would lead to an 
enhanced C3.4 coupling, the generation of a set of new op- 
erators to explain the polarization [|T], and at the same 
time not disturb all the successes of the standard model 
in the B and K sectors. Recent work on the polarization 
question was done in Ref . [p8| and |89j . 

Note: While this manuscript was in preparation 
Ref. |3^] appeared where Ab^ was also computed in the 
NDR scheme. We have verified that our result for Ab± 
in the NDR agrees with theirs as a function of u and z 
(they use a different basis, and the relevant comparison 
is for 2Ab [ P - 2Ac ( / ] ). Unlike Ref. ||, for the up and 
charm loops we demonstrated the simplicity of using the 
offshell UV subtraction procedure, and also presented the 



computation in the HV scheme for 75. Ref. |35J includes 

( f) • — 
small Cj>3a s terms in Ab\ , which we neglected because 

they are expected to compete with other terms of simi- 
lar numerical size (such as the complete two-loop correc- 
tions with Ci ) which remain unknown. We also derived a 
new factorization theorem for chiraly enhanced penguins, 
and demonstrated that only the short-distance perturba- 
tive coefficients give imaginary parts to the correspond- 
ing amplitudes. In contrast Ref. [j35| includes a complex 
hadronic parameter in their modeling of the analogous 
terms which they obtain with an additional expansion 
in a s (pn), where p 2 ~ EA is the intermediate scale. Fi- 
nally our phenomenological analysis differs from Ref. f35| . 
Our strategy was to avoid expanding in a s (pi), and to 
use data on the tree-amplitudes to determine the most 
important non-pcrturbative parameters. This allowed us 
to reduce the model uncertainty considerably while still 
predicting the penguin amplitudes. In contrast Ref. |35| ] 
models all non-perturbative parameters, and hence has 
larger parameter uncertainty in their final result. 
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APPENDIX A: ONE-LOOP FUNCTIONS 



In this appendix we quote some of the basic loop integrals used in the text. We work in d = 4 — 2e dimensions in 
the MS scheme, and the (47r) _1 e TE constant is absorbed in momentum subtraction scale p 2 . The basic loop integrals 
Jo and Io{q 2 ) are defined as, 



■h{p-,q) = 



16tt' 



d d k 



1 



2p-q 



Lio 



(2ir) d [k 2 - m 2 c ] [{k + p) 2 - m 2 ] [{k - q) 2 - m 2 ] 



Li 2 



1 - yjl - Am 2 Jq 2 
-Li 2 



16^ 2 



d d k 



1 + ^1 - 4m 2 c /(p + q) 2 
1 1 



1 + y/l - 4m 2 . /q 2 
Li 2 



1 - ^1 - 47712/(^ + ^)2 



(2n) d {k 2 - m 2 ) {{k-q) 2 ~m 2 c ) 



-+h(q 2 



2 + In (4) - % 2 -4m^)Vl- 47772/^2 / ln fi+Vl 



1 + v/1 - Am 2 c l<t 



y/l — Am 2 /q 2 
-2 9(4m 2 c -q 2 )^4ml/q 2 - 1 cot" 1 (^Am 2 c /q 2 - l) , 



(Al) 



Some limits for these functions are The combinations of these integrals that appear in the 

t\q 2 ) EE hm 7o(<Z 2 ) = - + 2 + ln(4)+^ ! 

h = Hm / ( g 2) = I +m (iL), 
q 2 ->o e \mj.J 

ifV) = 2 + ]n(£)+«r, 

2 

In = lnfArV (A2) 
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body of the paper are 



K{u,z,p) = m 2 c J (p,q) , 



Io(q 2 ) = G 



(A3) 



where u and z are momentum fractions of quarks in M\ 
and M2 respectively. Notice that functions h?(u,z) are 
finite and dimcnsionless. Change of the functional de- 
pendence from square of momenta to momentum frac- 
tions comes from the relations q 2 = rn?(l — u)z and 
(p + q) 2 = mf(l-u). 

We also list some useful Fierz relations. For conve- 
nience we define, 

P Xa = [qnliT a P L b v ] [d Tl 4P L q fl ] 
V Xa = [q n T a P L b v ] [dn^P R qu] . 

Fcirzing gives following formulas useful for simplification, 



e_^[df^P L b v ] 


[q n 7X^T a q n ] 


_ 1 (p\a 

3 V 


- v Xa ) 




[qnlt^qn] 


18 v 


- v Xa 


r bc [d n 4PLT c b t 


} [q n llT b q n ] 


= -I(P*M 


-V Xa ) 


[dniPLlillbv] 


[q n l±^T a q ri ] 


3 ' 




[dniPLlllib v ] 


[q n J±vT a q fl ] 


_ ^ pXa 

3 ' 




\d n 4,P L b % 


] [q nl iT a q fl ] 


= -\{P Xa l 


- v Xa ) 



(A4) 

These results include the minus sign from permuting 
fermion fields. 



APPENDIX B: CHIRALY ENHANCED TERMS 

IN SCETn 



In section VIII we made the statement that the op- 



erators Q\ in Eq. ([77]) only contribute for B* de- 
cays. The simplest way to verify this statement is to 
consider the Dirac structures generated by matching the 
T2 time-ordered product of these operators onto oper- 
ators in SCETn- This can be done working to all or- 
ders in a s . From Ref. || the most general perturbative 
matching generates Wilson coefficients given by jet func- 
tions J and Jj_ whose form is constrained by RPI, chiral- 
ity, power counting and dimensional analysis [u)\ = zuj, 



1 f . f dk + +ik + y -/ 2 
2tt 



= (l-z)u, X = l-X, Xn,uj = {Wl£n)<j], 

T [(ZnW) Ul igB^ 4 P R , L ] ia (0) [igfaW^ n ] j \y) 

= i S ab S(y+)S^(y ± )- I dx 
u Jq 

x{-J ± (z,x,k + ) (^P R ,Lll7l ) . . [x^HXn]-™] 

+ J(z, x, k+)(iP L , Rl l),pCn^ix^%J\}, (Bl) 

where {i,j} and {a, b} are spin and color indices. This re- 
sult on the RHS includes the sign from antipcrmuting the 
fermion fields that are contracted with the spin indices 
ij. To use this formula other Dirac structures occurring 
in the T2 time-ordered product should be grouped with 
the heavy-quark field h v and soft-quark field q s . Using 
Eq. (Bl) we find that all the operators Q^ 2x> * give the 
structure q s iji^ c ^h v which has a vanishing i?-mcson ma- 
trix element, but would be nonzero for B* initial states. 

In section |VIII we used C^ M ( Z ) ~ z which followed 
from the power counting in SCETi that indicates that it 
can not be more singular than (f M (z). We also stated 
that this scaling could be checked by factorizing the 
SCETi time-ordered product that defines this form fac- 
tor using SCETn- To do so we again use Eq. (Bl) and 
take the matrix element of the resulting operators to find 



/s/m 



Jo i z 

(B2) 



Here J±(z, x, k+) = S(x — z)na s (fi)CF/(N c xk + ) at low- 
est order, so the behavior of C^ M ( Z ) as z ^ is inher- 
ited from 4>pp(x) as x — > 0, giving linear scaling oc z. 
The one-loop result for J± is also known [il| and 
the linear scaling is reproduced at this order. The limit 
x — > corresponds to the collincar quark in the form fac- 
tor becoming soft, and for this matrix element there is 
no corresponding diagram with a soft quark at this order. 
Hence, because we do not expect an overlap with a soft 
diagram, we do not expect there is a need for any zero-bin 
subtractions, and hence no endpoint divergences which 
would result from constant scaling as z — > 0. The other 
interesting limit is z — ► 1, where we expect C^ M ( Z ) ~ 1- 
From the point of view of SCETn this limit is more in- 
teresting because there are diagrams with soft antiquarks 
in the form factor, and we must avoid double counting 
them. Indeed the tree level jet function appears to give 
C^ M (z) ~ 0pp(z)/(l — z) 2 , which would imply singular 
behavior as z — > 1. However in SCETn to avoid dou- 
ble counting the region where this quark is soft we must 
make zerobin subtractions in defining this singular mo- 
ment |f30|| . These subtractions modify the distribution, 
causing scaling behavior of <f)pp bm (z) ~ (1 — z) 2 in the 
endpoint region (and also generate dependence of this 
distribution on an additional rapidity parameter). The 
result is that as z — > 1, (z) ~ 1 as expected from 
power counting in SCETi. 
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APPENDIX C: BETA'S AT (i = m b /2 AND /i = 2m b also required the values at /i = m fc /2 = 2.3 GeV 



In Eq. (119) we quoted values for the annihilation mo- 
ments at /i = nib = 4.7 GeV. For our error analysis we 



$™ = (-5.5 ± 3.0) x 1(T 2 , 

= -2.35 ±0.78, 
= 0.0 ±9.7, 

/3jf = -0.27 ±0.11, 

PI* = -2.45 ± 0.80 , 

/3£f = 0.0 ± 12.3, 

/3f c ^(9.2^)xl0- 3 , 

Ci = (-7-0^i)xl0- 2 , 

and at /i = 2rrib 

= (-1.6 ±0.9) x 10~ 2 , 
/9j& = -0.73 ±0.23, 

= 0.0 ±2.7, 
PI* = -0.095 ± 0.069 , 
/9jf£ = -0.76 ±0.24, 
/3jf = 0.0 ± 3.4 , 
^ = (2.7t 2 ; 2 7 )xl0- 3 , 
C 1 = (-2.2ll;»)xl0- a > 



= 1.15±0.58, 
022, =0.61 ±0.34, 
=0.0 ±9.1, 

^ = 0.61 ±0.34, 

= 0.0 ± 11.2, 
^ = -0.19^, 
C 2 = (-1.411;?) xlO" 2 , 

= 0.35±0.18, 
/3hc2 = -0.015 ± 0.040 , 
= 0.0 ± 2.5 , 

^1 = -0.015 ± 0.040 , 
fig = 0.0 ± 3.0 , 
/3£ = -0.06t;°2, 
C 2 = (L7+i3) xlO" 3 , 



= -0.25±0.15, 
/3^3 = (4-4±5.3) x 10" 3 , 
= 0.0 ±0.13, 

/?Ll = (0.8±5.5)xlO- 3 , 
/2jf = 0.0 ±0.17, 
/4 P =(4.3i 3 j)xl0- 2 , 
C3 = (-l-4±r: 2 o)xl0- 4 , 

Z?™ = -0.089±0.054, 
/?L 7r 3 = (-4.9±0.8)x lO- 3 , 
0$ = 0.0 ± 0.040 , 

^ = (-4.9±0.8)x 10" 3 , 
/3jf = 0.0 ± 0.051 , 

/3 4 P c P = (l-5t^)xlO- 2 , 
C 3 = (-2-l±l-7)xl0- 4 , 



(CI) 

/3S = (-8.6±2.8)x 10~ 2 , 
/3££ = 0.0 ±0.16, 

^5 = (-8-6 ± 2.7) x 10- 2 , 
= 0.0 ± 0.20 , 

^4-(2.6t^)xl0- 3 , 



(C2) 

/3^ 4 = (-3.7±l.l)xl0" 2 , 
(3^ = 0.0 ± 0.053 , 

^ C 4 = (-3.7±l.l)xl0" 2 , 
/3jf = 0.0 ± 0.065 , 

C 4 = (l.H^)xl0- 3 . 
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